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Annual  Technical  Report 


Introduction 


In  this  report,  research  In  seismology  on  fundamental  prob¬ 
lems  related  to  the  detection,  location,  and  Identification  of 
underground  nuclear  explosions  Is  described.  The  period  covered 
by  the  report  is  September,  1967  through  August,  1968. 

The  personnel  engaged  In  this  work  are  all  associated  with 

t 

the  Department  of  Geophysics  and  Geophysical  Engineering,  Instl- 
t  tute  of  Technology,  Saint  Louis  University.  Scientists  In  charge 

of  work  on  various  tasks  under  the  project  are:  S.  J.  Duda,  C.  Kiss- 
llnger,  0.  W.  Nuttll,  and  W.  V.  Stauder.  In  addition,  ten  research 
assistants,  all  graduate  students  In  geophysics,  contributed  to 
the  research. 

In  the  report,  the  completion  of  a  study  of  Love  and  Rayleigh 
wave  propagation  In  the  central  part  of  North  America  Is  described. 
Observational  and  theoretical  studies  of  seismic  body  waves,  P  and 
S,  follow,  and  finally,  research  In  model  seismology  Is  described. 

A.  Work  in  Progress. 

1.  Long  Period  P-Wave  Amplitudes  at  LA3A . 

The  object  of  this  study  Is  to  Investigate  the  amplitude 
anomalies  of  the  long  period  P  waves  at  LASA  as  a  function  of  dis¬ 
tance  and  azimuth  to  the  source  region,  with  emphasis  upon  the 
Interpretation  of  the  anomalies  in  terms  of  crust  and  upper  mantle 
structure . 
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Two  approaches  to  this  problem  have  been  taken.  One  compares 
P  amplitudes  at  the  individual  stations  for  a  given  earthquake,  to 
determine:  l)  if  any  particular  station  has  consistently  larger 
or  smaller  amplitudes  relative  to  the  other  stations;  2)  if  there 
is  an  amplitude  variation  across  the  array  which  is  a  function  cf 
distance  and  azimuth;  and  3)  If  there  are  any  epicentral  distances 
for  which  the  change  of  amplitude  with  distance  i3  abnormally 
large,  indicative  of  rapid  velocity  change  with  depth  in  the  mantle. 
The  second  approach  uses  the  departure  of  the  P  motion  from  the 
great  circle  path  as  determined  from  the  amplitudes  of  the  hori¬ 
zontal  components  of  ground  motion  at  the  individual  stations  of 
the  array.  These  departures  from  the  great  circle  path,  or  azi¬ 
muth  anomalies,  are  theoretically  related  in  a  relatively  straight¬ 
forward  manner  to  the  lateral  changes  in  velocity  that  may  exist 
beneath  LASA. 

Much  effort  was  devoted  to  the  latter  approach,  particularly 
in  the  gathering  of  data,  because  the  method  has  the  potential  of 
yielding  in  a  direct  manner  information  about  the  dip  or  horizon¬ 
tal  velocity  changes  of  the  subsurface  layers.  The  data,  however, 
exhibited  an  amount  of  scatter  equal  at  least  in  size  to  any 
existing  anomaly,  and  attempts  to  reduce  the  scatter  proved  fruit¬ 
less.  For  this  reason  that  part  of  the  study  has  been  dropped. 

P  wever,  because  no  similar  studies  have  been  reported  in  the 
literature,  some  discussion  of  the  difficulties  encountered  is 
included  here. 

The  basic  data  were  copies  of  the  Develecorder  films  of  th° 
long  period  seismograms  of  the  D  and  F  rings.  Earthquakes  which 
produced  impulsive  P  waves  at  LASA,  xvere  selected  with  the 


additional  requirements  that  the  mlcroselsinlc  background  level 
was  low  and  that  the  P  amplitudes  were  large  but  not  so  large 
as  to  have  the  turning  points  of  the  motion  fall  off  the  film. 

In  the  first  attempt  to  obtain  the  necessary  data  the  amplitudes 
of  the  zero-to-flrst  maximum  (or  minimum)  and  the  first  maximum- 
to-first  minimum  were  read  with  the  assistance  of  a  film  digitizer. 
These  horizontal  amplitude  components  were  combined  to  give  an 
observed  direction  of  wave  arrival  (0)  and  compared  with  the  com¬ 
puted  azimuth  (C)  of  the  great  circle  path  at  the  station.  Re¬ 
sults  for  some  of  the  earthquakes  studied  are  given  in  Table  1. 


Table  1 

Earth¬ 
quake  Origin  Sta-  0  *  0  **  0,-C  0o-C 

Date  Time  Lat.  Long,  tion  1  2  C  1  2 

13  May  67  05-18-55.4  56. 5N  152. 6W  D-l  314.1°  309.3°306.1°  +8.0°  +3.2° 

D-2  315.0  308.8  306.6  +8.4  +2.2 

D-3  311.2  301.4  306.5  +4.7  -5.1 

P-3  —  309.6  307.3  —  +2.3 

P-4  --  310.5  305.3  —  +5.2 

21  May  67  07-18-12.8  27. 9N  111.3W  D-l  195.3  194.8  194.6  +0.7  +0.2 

D-2  202.4  201.4  194.5  +7.9  +6.9 

D-3  194.9  187.9  193.2  +1.7  -5.3 

D-4  195.8  —  193.2  +2.6 

P-1  197.3  199.2  196.0  +1.3  +3.2 


P-2  191.4 

203.5 

196.4 

-5.0 

+7.1 

P-3  193.2 
P-4  194.2 

191.0  191.9 

+1.3 

-0.9 

01  Jul  67  23-10-07.2 

54. 4N  158. OW 

192.2 

191.5 

+2.7 

+0.7 

D-l  296.1 

-- 

302.7 

-6.6 

-  - 

D-2  298.4 

— 

303.1 

-4.7 

D-3  292.4 

-- 

302.9 

-10.5 

— 

D-4  296.2 

— 

302.4 

-  6.2 

— 

29  Jul  67  10-24-24.6 

P-4  297.2 

— 

301.7 

-4.5 

— 

6.8N  73. OW 

D-l  137.1 

-- 

134.3 

+2.8 

— 

D-2  133.3 

-- 

133.9 

-0.6 

-- 

P-3  130.3 

132.3 

-2.0 

-- 

P-4  137.3 

w 

133.4 

+3.9  ' 

*  Prom  first  half-amplitude 


**  From  first  peak  to  peak  motion 
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No  pattern  of  regularity  appeared  in  the  0-C  data  of  Table  1 
and  in  other  similar  data  not  included  in  this  report.  In  fact, 
when  these  0-C  data  were  plotted  against  the  azimuth  of  the  great 
circle  path  it  was  found  for  each  of  the  azimuths  for  which  there 
were  sufficient  data  (130°,  190°,  and  310°)  that  the  negative  and 
positive  values  of  the  0-C  were  about  equal  in  value  and  in  num¬ 
ber.  This  indicates  that  the  cause  of  the  anomalies  is  some  type 
of  reading  or  instrumental  effect  rather  than  lateral  velocity 
variations  in  the  crust  and  mantle. 

Although  the  earthquakes  selected  for  study  had  relatively 
small  microseismic  and  instrumental  noise  amplitudes,  it  is  pos¬ 
sible  that  a  peak  in  noise  motion  could  arrive  simultaneously 
and  interfere  with  a  peak  in  the  P  wave  motion.  Azimuths  calcu¬ 
lated  from  peak  amplitudes  of  the  P  motion  would  thus  be  in  error. 
This  difficulty  can  be  obviated,  however,  by  the  use  of  particle 
motion  diagrams  if  the  amplitude  of  the  noise  is  small  compared 
to  that  of  the  signal  and  if  both  do  not  have  the  same  frequency 
content . 

Figure  1  contains  horizontal  component  particle  motion  dia¬ 
grams  for  the  P  waves  of  21  May  1967  as  recorded  at  the  D-l  and 
F-4  stations.  It  can  be  seen  from  this  figure  that  the  motion 
is  approximately,  but  not  exactly,  linear.  As  a  result  of  this 
departure  from  linearity  there  is  uncertainty  in  the  value  assigned 
to  the  direction  of  the  particle  motion.  Table  2  contains  the 
direction  of  a  visual  estimate  of  a  straight  line  fit  to  the  mo¬ 
tion  in  the  first  and  third  quadrants,  and  the  corresponding 
departure  from  the  great  circle  path  for  all  the  stations  of 
the  D  and  F  rings  for  the  earthquake  of  21  May  1967. 
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Station 

V 

Table  2 

°rc 

0  ** 
III 

°III 

D-l 

197.4° 

+2.8° 

196.9° 

+2.3 

D-2 

198.4 

+3.9 

200.9 

+6.4 

D-3 

191.9 

-1.3 

193.1 

-0.1 

D-4 

190.5 

-2.7 

188.4 

-4.8 

F-l 

190.8 

-5.2 

194.4 

-1.6 

F-2 

204.5 

+8.1 

202.1 

+5.7 

P-3 

187.4 

-4.5 

193.2 

+1.3 

F-4 

188.5 

-3.0 

193.8 

+2.3 

*  Prom  motion  in  first  quadrant. 

**  Prom  motion  in  third  quadrant. 

Finally,  to  remove  the  subjectivity  involved  in  assigning  a 
direction  to  the  particle  motion  "by  eye, "  a  least-squares  fit  of 
a  straight  line  was  made  to  the  data .  The  results  are  summarized 
in  Table  3. 

Table  3 


Station 

V 

°rc 

°2** 

°2'C 

D-l 

189.4° 

-5.2° 

195.6° 

+1.0 

D-2 

197.9 

+3.4 

196.5 

+2.0 

D-3 

196.7 

+3.5 

193.4 

+0.2 

D-4 

191.1 

-2.1 

187.9 

-5.3 

F-l 

198.7 

+2.7 

192.4 

-3.6 

F-2 

194.9 

-1.5 

203.6 

+7.2 

P-3 

198.2 

+6.3 

189.2 

-2.7 

F-4 

196.8 

+5.3 

190.8 

-0.7 

♦From  first  half-amplitude. 

♦♦Prom  first  peak  to  peak  motion. 

The  values  contained  in  Tables  1  and  3  should  be  similar, 
although  comparison  indicates  no  such  agreement,  even  insofar  as 
the  sign  of  the  O-C's  is  concerned.  Because  the  values  given  in 
Table  2  were  determined  from  several  cycles  of  the  P  motion,  they 
are  expected  to  be  somewhat  different  than  those  in  Tables  1  and 
3.  In  any  case,  if  dip  of  the  M  discontinuity  or  lateral  velocity 
variations  in  the  mantle  under  LASA  were  responsible  for  the 
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azimuth  anomalies,  then  for  a  given  earthquake  rll  the  O-C's 
should  have  the  same  sign  and  have  about  the  same  value  at  all 
the  LASA  stations.  It  can  be  seen  from  the  tables  that  positive 
and  negative  values  of  0-C  for  the  earthquake  of  21  May  1967  are 
about  equal  In  number. 

It  Is  concluded  that  the  observed  azimuth  anomalies  are 
most  probably  caused  by  some  combination  of  mlcroselsms.  Instru¬ 
mental  noise,  local  and  very  shallow  station  site  Irregularities, 
and  possibly  small  errors  In  Instrument  calibration  and  seis¬ 
mometer  orientation.  Their  combined  effects  appear  to  produce 
an  apparent  azimuth  anomaly  of  a i  much  as  5°  In  a;  solute  value. 

The  second  approach  to  amplitude  anomaly  studies  of  long 
period  P  at  LASA  has  proved  more  fruitful,  even  though  the  work 
Is  in  a  preliminary  state.  One  Important  conclusion  Is  that  the 
long  period  P  wave  closely  maintains  Its  waveform  and  amplitude 
across  LASA,  with  departures  of  Individual  station  amplitudes 
from  the  mean  amplitude  usually  no  greater  than  ten  percent. 

This  suggests  that  long  period  P  waves  are  Ideally  suited  for 
such  processing  techniques  as  delay  and  sum,  etc.,  whenever  the 
seismic  source  Is  of  sufficiently  large  magnitude  to  generate 
them. 

Table  4  contains  the  hypocentral  coordinates  of  nine  earth¬ 
quakes  used  In  this  part  of  the  study.  The  eplcentral  distance, 
azimuth,  and  vertical  component  P  wave  amplitude,  as  well  as  the 
departure  from  the  average  amplitude  of  an  Individual  station's 
amplitude,  are  given  for  the  nine  earthquakes  In  Table  5.  The 
amplitudes  given  In  Table  5  are  peak  to  peak  values,  measured 
from  the  first  maximum  to  the  first  minimum,  or  vice  versa. 
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Table  4 


Earth- 


quake 

Number 

Date 

Origin 

Time 

Lat . 

Long. 

Depth 

Magnl 

tude 

1 

21  May  67 

07-18-12.8 

27. 9N 

111.3  V/ 

33  ton 

4.7 

2 

18  Oct  67 

01-11-44.8 

79. 8N 

2.4E 

33 

33 

5.7 

3 

13  May  67 

05-18-55.4 

56. 5N 

152.61/ 

5.3 

4 

27  May  67 

17-22-58.7 

51. 9N 

176. IE 

34 

5.8 

3 

16  May  67 

12-58-09.5 

13. 5N 

90.6V/ 

95 

4.8 

6 

19  Jun  67 

17-07-45.4 

52. 7N 

166.9W 

33 

5.7 

7 

10  Dec  67 

12-06-50. 3 

40. 5N 

124. 6W 

5 

5.8 

8 

20  May  68 

20-05-49.06 

30.7S 

178.41/ 

46 

6.0 

9 

21  May  68 

08-20-00.86 

44. 9N 

150. IE 

noi  mal 

5.7 

Table  5 


Earthquake 

Z  amplitude 

departure 

Station 

Number 

Distance 

Azimuth 

(milimicrons) 

from  mean 

D-l 

1 

19.4° 

194.6° 

37.5 

-10.7 

2 

47.4 

13.3 

91.8 

-  7.3 

3 

30.0 

306.1 

13.2 

-20.9 

4 

48.8 

306.3 

49.9 

-  6.9 

5 

35.6 

153.9 

15.1 

+12.7 

6 

38.8 

301.8 

23.1 

-11.1 

7 

14.9 

251 .7 

107.5 

-14.5 

8 

101.0 

236.8 

172.8  (?) 

-29.6  (?) 

9 

66.8 

311.3 

35.1 

-  1.9 

D-2 

1 

19.0 

194.5 

39.3 

-  6.4 

2 

47.7 

13.2 

94.3 

-  4.7 

3 

30.1 

306.6 

14.6 

-12.6 

4 

48.9 

306.5 

51.4 

-  4.1 

5 

35.4 

153.5 

12.2 

-  8.9 

6 

38.9 

302.1 

23.6 

-  9.2 

7 

14.7 

252.7 

— 

8 

100.8 

236.7 

271 .8 

+10.7 

9 

66.9 

311.4 

35.6 

-  0.5 

D-3 

1 

19.0 

193.2 

42.1 

+  0.2 

2 

47.8 

13.2 

105.9 

+  1.5 

3 

29.8 

305.5 

18.0 

+  7.8 

4 

48.6 

306.3 

57.0 

+  6.3 

5 

35.6 

152.7 

10.2 

-23.9  (?) 

6 

38.6 

302.0 

28.3 

+  8.8 

7 

14.4 

251.8 

117.4 

-  6.6 

8 

100.5 

236.4 

-- 

9 

66.7 

311.1 

-- 

—  — 

D-4 

1 

19.4 

193.2 

44.2 

+  5.2 

2 

47.4 

13.3 

99. 5 

+  0.5 

3 

29.6 

305.0 

19.6 

+17.4 

4 

48.4 

306.0 

57.8 

+  7.8 

5 

35.9 

153.1 

13.7 

+  2.2 

6 

38.4 

301.6 

27.1 

+  4.2 

7 

14.6 

250.5 

116.4 

-  7.4 

8 

100.8 

236.5 

241.7 

-  1.5 

9 

66.5 

311.0 

^5-8 

0.0 

iftw— 
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Earthquake 

Z  amplitude 

%  departure 

Station 

Number 

Distance 

Azimuth 

(milimicrons) 

from  mean 

P-1 

1 

20.0 

196.0 

42.4 

+  0.9 

2 

46.7 

13.5 

— 

3 

30.1 

305.4 

14 .6 

-12.6 

4 

48.8 

306.2 

59.9 

+11.7 

5 

35.9 

155.3 

14.6 

+  8.9 

6 

38.9 

301.4 

31.9  4 

+22.7 

7 

15.5 

250.9 

168.1  (?) 

+33.7  ( 

8 

101.7 

237.4 

255.5 

+  4.1 

9 

66.8 

3H.5 

39.5 

+10.3 

P-2 

1 

18.5 

196.4 

39.8 

-  5.2 

2 

48.2 

13.1 

92.7 

-  6.4 

3 

30.8 

307.4 

15.8 

-  5.4 

4 

49.5 

307.1 

48.7 

-  9.1 

3 

3^.7 

153.9 

13.5 

+  0.7 

6 

39.5 

302.9 

19.1 

-26.5 

7 

14.9 

255.6 

134.5 

+  7.0 

8 

100.8 

236.9 

229.1 

-  6.6 

9 

67.6 

311.8 

29.5 

-17.6 

P-3 

1 

18.3 

191.9 

44.3 

+  5.4 

2 

48.4 

13.0 

98.8 

-  0.2 

3 

29.8 

307.3 

18.0 

+  7.8 

4 

48.6 

305.5 

52.9 

-  1.3 

5 

35.2 

151.4 

13.9 

+  3.7 

6 

38.5 

302.4 

23.8 

-  8.5 

7 

13.9 

253.0 

130.6 

+  3.9 

8 

99.9 

235.9 

252.6 

+  2.9 

9 

66.7 

311.0 

34.8 

-  2.8 

P-4 

1 

19.8 

191.5 

46.3 

+10.2 

2 

47.0 

13.3 

106.4 

+  7.5 

3 

29.1 

305.3 

20.2 

+20.9 

4 

47.8 

305.5 

51.4 

-  4.1 

5 

36.5 

152.6 

13.7 

+  2.2 

6 

37.9 

300.9 

31.9 

+22.7 

7 

14.4 

247.9 

105.8 

-15.8 

8 

100.8 

236.2 

295.1  (?) 

+20.2  ( 

9 

65.9 

310.6 

40.4 

+12.8 

Prom  Table  5  It  can  be  observed  that  stations  D-l,  D-2,  and 
P-2  exhibit  a  tendency  to  have  amplitudes  less  than  the  average, 
whereas  D-3,  D-4,  P-1,  and  P-4  have  amplitudes  greater  than  the 
average.  These  observations,  however,  are  based  upon  a  small 
amount  of  data  and  may  be  changed  as  more  data  accumulate. 

There  are  at  present  insufficient  data  to  determine  if  the 
amplitudes  have  any  azimuthal  dependence.  For  the  epicentral  dis¬ 
tance  intervals  studied  there  is  no  indication  of  a  rapid  decrease 
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or  Increase  In  amplitude  with  Increasing  eplcentral  distance, such 
as  are  produced  by  discontinuous  Increases  In  the  velocity  gradi¬ 
ent  In  the  mantle.  However,  as  yet  data  Is  lacking  for  some  epl¬ 
central  distances. 

To  repeat,  the  most  Important  conclusion  to  be  drawn  from 
the  amplitude  observations  to  date  Is  the  coherence  of  the  long 
period  P  waveform  across  LASA,  with  the  amplitude  of  P  at  an  in¬ 
dividual  station  usually  differing  by  no  more  than  ten  percent 
from  the  average  of  the  values  at  all  the  stations  of  the  D  and 
F  rings. 

2.  Body  Wave  Travel  Time  Anomalies 
Statement  of  the  Problem. 

Travel-time  tables  cf  body  waves  (e.g.  Jeffreys-Bullen)  bear 
Information  concerning  the  average  structure  of  the  Earth.  Travel 
time  anomalies,  observed  at  various  places  over  the  Earth,  indi¬ 
cate  heterogeneities  in  the  Earth’s  structure.  The  problem  in¬ 
vestigated  is  1)  to  measure  the  travel  times  and  find  statistically 
significant  deviations  from  the  mean,  2)  to  determine  the  probable 
location  of  the  anomalous  structure,  and  3)  to  interpret  the  anom¬ 
alies  from  a  tectonophyslcal  point  of  view. 

In  return,  it  is  expected  that  the  more  detailed  knowledge 
of  the  travel-time  anomalies  and  the  structure  of  the  Earth  v/ill 
increase  the  accuracy  of  epicenter  determination  of  natural  and 
artificial  events  on  the  Earth. 

In  continuation  of  previous  efforts  (see  Annual  Technical 
Report  Contract  AF  19(628) -5100,  September  1967,  (ATR  1967),  P 
Time  Delays  and  Continental  Uplift)  the  investigation  has  been 
extended  in  the  following  directions: 
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1.  determination  of  the  first  derivatives  of  the  travel-time 
curve  from  observations  of  P-arrlvals  In  the  area  pre¬ 
viously  found  to  have  anomalous  relative  P-tlme  delays 
(RPP); 

2.  finding  of  a  theoretical  model  of  crust  and  upper  mantle 
yielding  relative  P-tlme  delays  and  first  derivatives  of 
the  travel-time  curve  consistent  with  the  observed  ones; 

3.  application  of  the  previous  approaches  with  some  modifica¬ 
tions  to  S-v/ave  travel  times. 

Determination  of  first  derivatives  of  P-wave  travel  times. 

Prom  the  arrival  times  of  P-waves  at  two  stations  and  the 
epicentral  distances  of  the  earthquake,  the  derivative  of  the 
travel-time  curve  can  be  computed.  If  regional  anomalies  are 
being  studied,  the  stations  must  not  only  have  similar  epicentral 
distances,  but  also  be  situated  close  to  each  other.  The  deriva¬ 
tive  then  characterizes  the  structure  between  the  two  stations. 

In  Pennoscandia,  the  region  selected  for  the  research  for  reasons 
outlined  in  ATR  1967,  inhomogeneities  in  the  structure  seem  to 
exist  over  relatively  short  distances .  If  the  distance  between 
the  two  stations  is  too  large,  the  anomaly  in  the  derivative, 
if  existent,  would  lose  its  meaning.  Therefore,  the  derivatives 
were  not  computed  for  all  possible  combinations  of  15  stations 
into  pairs,  but  were  computed  for  all  station  pairs  for  which 
the  distance  between  the  stations  did  not  exceed  about  300  ml. 

The  derivatives  were  classified  according  to  the  azimuths  toward 
epicenters  into  four  groups  (N,  E,  S,  W),  and  displayed  as  func¬ 
tions  of  epicentral  distance.  Due  to  the  natural  distribution 
of  sei8mically  active  regions,  most  earthquakes  arrived  from 
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northerly  and  easterly  directions,  and  the  derivatives  in  those 
groups  (N,  E)  are  of  primary  interest.  The  derivatives  in  the 
other  two  groups  (S,  Vi)  are  not  used  as  independent  evidence, 
but  only  as  supporting  the  conclusions.  Prom  a  detailed  analysis 
of  all  the  derivatives  it  was  found  that  for  some  station  pairs 
the  derivative  deviates  significantly  either  in  the  positive  or 
negative  direction,  while  for  some  station  pairs  no  deviation  is 
found.  Pig.  2  gives  samples  of  all  three  cases* 

No  deviation  of  the  derivative  from  the  expected  value 
means  that  there  is  no  lateral  change  in  the  structure  between 
the  station  pair,  but  an  anomalous  structure)  unifobm  over  the 
distance  between  the  stations,  cannot  be  excluded «  A  positive 
or  negative  deviation  of  the  derivative  means  a  lateral  change 
in  the  structure  between  the  stations.  Here  it  is  implied  that 
the  source  of  the  anomaly  in  the  derivative  of  the  travel-time 
curve,  as  for  the  RPD’s  (see  ATR  1967) »  lies  in  the  uppermost 
part  of  the  earth,  beneath  the  station.  This  assumption  seems 
justified  for  at  least  two  reasons: 

i.)  Measurable  lateral  inhomogeneities  in  the  Earth  are  more 
probable  at  shallower  depths  than  at  greater  depths  be¬ 
cause  of  the  increasing  role  of  hydrostatic  pressure  at 
greater  depths,  tending  to  create  spherically  symmetric 
conditions  around  the  center  of  the  Earth,  and  because 
of  tectonic  forces  and  tectonic  structures  known  with 
some  certainty  to  be  inhomogeneous  down  at  least  to  the 
Moho-discontinuity,  and  probably  deeper, 
il.)  The  two  seismic  rays  originating  at  the  earthquake  focus 
and  terminating  at  the  station  pair  under  consideration 
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diverge  from  each  other  mainly  just  prior  to  the  arrival 
at  the  stations.  Thus,  beneath  the  stations  they  are 
most  likely  to  traverse  anomalous  structures  and  suffer 
a  relative  acceleration  or  deceleration. 

This  assumption  Is  supported  by  the  evidence  of  Pig.  2.  SKA, 

UME,  and  KJN  are  three  of  the  studied  stations,  situated  on  a 
roughly  E-Vi  striking  line  crossing  Pennoscandla .  The  epicenters 
in  group  N  lie  mainly  in  the  Aleutian  Islands,  north-northeast 
of  Pennoscandla.  The  derivatives  between  SKA-UME  and  KJN-UME 
were  computed.  The  arrival  times  produce  a  distinct  positive 
anomaly  at  SKA-UME  and  a  distinct  negative  one  at  KJN-UME.  Sim¬ 
ilarly,  the  epicenters  in  group  E  produce  a  slightly  positive 
anomaly  at  SKA-UME  and  a  clear  negative  one  at  KJN-UME.  It  is 
also  interesting  to  note  that  in  SKA-UME,  while  the  anomaly  for 
earthquakes  from  easterly  directions  (group  E)  are  positive,  as 
mentioned,  the  anomaly  for  earthquakes  from  westerly  directions 
(group  17)  are  negative,  as  far  as  can  be  seen  from  the  limited 
number  of  epicenters  in  the  last  group.  These  features  will  be 
explained  below  with  the  help  of  a  postulated  model. 

Theoretical  model  of  crust  and  upper  mantle. 

The  derivatives  of  the  travel -time  curve  for  various  epi- 
central  distances  X'jere  computed  and  analyzed  for  a  total  of  20 
station  pairs.  The  findings  are  concordant  with  those  of  an 
analysis  of  the  relative  P-time  delays  insofar  as  higher  veloc¬ 
ities  in  the  central  portion  of  Pennoscandla  are  indicated.  In 
addition,  the  derivatives,  because  of  closer  spacing  of  the  obser¬ 
vation  points,  permit  the  delineation  of  the  limits  of  regions 
of  anomalous  velocities.  A  corresponding  contour  map  is  under 


construction. 


As  the  basis  of  the  analysis  a  layer  of  constant  thickness 
with  seismic  velocity  Increasing  both  with  depth  and  horizon¬ 
tally  Is  considered.  A  vertical  cross-section  parallel  to  the 
horizontal  velocity  gradient  Is  shown  In  Pig.  3.  Two  stations 
are  thus  situated  at  points  vrlth  different  velocities.  The  ver¬ 
tical  velocity  Increase  Is  the  same  over  the  entire  layer  and 
will  not  cause  any  RPD's  different  from  zero  or  anomalies  In  the 
derivatives  of  the  travel-time  curves.  The  velocity  In  the  ver¬ 
tical  direction  Is  then  assumed  constant  and  equal  to  an  average 
value.  The  velocity  change  In  the  horizontal  direction  Is  re¬ 
sponsible  for  the  RPD's  different  from  zero  and  for  anomalous 
derivatives.  The  effect  of  a  linear  horizontal  velocity  Increase 
upon  the  derivatives  was  computed  for  different  spacings,  DX, 
between  and  positions,  XI,  of  the  stations,  for  different  thick¬ 
nesses  H  and  different  rates  of  Increase  "a"  of  the  velocity  In 
the  horizontal  direction.  A  sample  computation  Is  presented  In 
Pig.  4,  In  which  the  anomalies  of  the  first  derivative  of  the 
P-travel-tlme  curve  as  a  function  of  eplcentral  distance  under 
varying  specified  conditions  are  shown.  A  comparison  of  the  ob¬ 
served  and  computed  anomalies  Indicates  that  the  thickness  of 
the  layer  with  a  horizontal  velocity  Inhomogeneity  is  relatively 
thin,  of  the  order  of  50  km.  Under  this  circumstance  RPD's  re¬ 
sult  with  maximum  values  and  variations  with  azimuth  are  In  agree¬ 
ment  with  earlier  observations.  Accordingly,  it  is  considered 
most  likely  that  a  horizontal  velocity  gradient  exists  beneath 
Scandinavia . 

Consideration  of  possible  causes  of  the  anomaly  has  led  to 
the  conclusion  that  it  may  be  related  to  the  geologically  recent 
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glaciations  of  Pennoscandla .  The  extent  of  the  anomaly  suggests 
that  the  portions  most  heavily  glaciated  suffered  a  compression 
of  the  strata.  The  removal  of  the  ice  resulted  in  an  elastic, 
instantaneous  rebound  and  an  imperfect  elastic  rebound  not  com¬ 
pleted  yet,  as  revealed  by  the  present  continental  uplift  of 
Pennoscandla . 

This  interpretation  of  the  continental  uplift  as  time- 
dependent  dilatation  is  even  more  plausible,  as  newer  gravi¬ 
metric  observations  in  Pennoscandla  are  not  in  support  of  the 
once  much  favored  interpretation  of  the  uplift  as  isostatlc  re¬ 
adjustment.  A  recently  published  map  of  isostatlc  gravity  anom¬ 
alies  for  Sweden  does  not  indicate  any  mass  deficit  under  the 
country,  as  required  for  the  latter  interpretation.  On  the 
other  hand,  the  interpretation  of  the  uplift  as  dilatation  is 
in  agreement  with  these  seismologlcal  observations  as  much  as  it 
is  with  known  gravimetric,  geological  and  mareographic  data. 
S-tlme  delays. 

This  travel-time  study  is  of  special  interest  in  connection 
with  the  establishment  of  a  seismic  array  station  in  Norway  un¬ 
der  Project  VELA-Uniform.  The  work  on  P-times  is  essentially 
completed  and  it  will  be  supplemented  with  a  corresponding  study 
of  the  S-times.  Some  work  has  been  performed  already. 

The  P-times  used  in  the  study  were  reported  in  the  bulle¬ 
tins  usually  with  an  accuracy  to  0.1  sec.  The  S-times  are  usu¬ 
ally  reported  with  an  accuracy  of  only  1.0  sec.  Moreover,  the 
identification  of  the  S-phase  is  more  difficult,  which  leads 
eventually  to  systematic  errors  in  the  S-times.  Prom  this  it 
is  evident  that  in  order  to  obtain  statistically  significant 
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data  a  larger  amount  of  data  must  be  used,  besides  the  elimina¬ 
tion  of  systematically  errant  observations. 

Arrival  times  of  S-waves  at  16  Scandinavian  stations  from 
more  than  1200  earthquakes  have  been  collected.  Several  computer 
programs  for  the  calculation  of  relative  S-time  delays,  slopes  of 
the  S-travel  time  curve  and  related  quantities  have  been  adapted 
or  prepared  and  test-run.  Work  on  this  problem  is  under  progress. 

3.  S  Wave  Studies 

All  studies  of  S  waves  carried  out  in  the  past  year  under 
this  contract  had  the  purpose  of  determining  the  crust  and  mantle 
velocity  distribution  and  its  regional  variations.  This  informa¬ 
tion  is  required  for  the  accurate  location  of  seismic  events  and 
the  determination  of  their  magnitude,  as  well  as  for  predicting 
the  travel  time  and  amplitudes  of  seismic  waves  which  result  from 
a  hypothetical  earthquake  or  explosion. 

The  S  wave  studies  may  be  divided  into  three  parts:  travel 
times  and  amplitudes  of  S  generated  by  nuclear  explosions  and 
earthquakes  as  recorded  by  seismograph  stations  throughout  North 
America,  apparent  wave  slownesses  (dT/dL  )  of  S  waves  across  LASA, 
and  studies  of  the  spectral  properties  of  the  S-wave  ground  mo¬ 
tion  and  their  relation  to  crustal  structure. 

Travel  Times  and  Amplitudes. 

During  the  twelve-month  period  covered  by  this  report  work 
was  completed  on  travel-time  and  amplitude  studies  of  S  for  a 
surface  focus  as  recorded  in  North  America  at  distances  of  1700 
to  6800  km.  This  work  complements  earlier  determinations  of  ob¬ 
servational  travel  times  for  the  same  distances  for  focal  depths 
of  33  and  120  km.  By  making  use  of  all  three  travel-time  curves 


we  have  calculated  an  S  wave  velocity  model  for  the  upper  1500  km 
of  the  mantle  which  Is  representative  of  an  "average  United  Stdtes" 
mantle  structure.  These  data  establish  the  existence  of  a  low- 
velocity  channel  at  a  depth  of  150  to  200  km,  a  discontinuous 
velocity  Increase  at  400  km  depth,  and  a  large  but  continuous 
velocity  Increase  at  750  km.  The  velocity  gradients  are  almost 
zero  In  the  depth  Intervals  50  to  150,  250  to  350,  and  500  to  70C 
km.  Some  of  the  data  appear  to  demand  In  addition  the  existence 
of  secondary  low- velocity  channels  at  depths  of  340  to  370  and 
670  to  710  km,  although  the  evidence  is  not  conclusive.  In  an 
attempt  to  resolve  this  important  question  data  which  will  make 
possible  the  construction  of  observational  travel-time  curves 
of  S,  sS  and  ScS  for  a  focal  depth  of  600  km  are  being  gathered. 
Data  extending  the  travel-time  curves  of  S  for  depths  of  33  and 
120  km  to  distances  beyond  10,000  km  for  the  purpose  of  obtaining 
lower  mantle  velocity  information  are  also  being  accumulated. 

Figure  5  compares  the  travel-time  curve  developed  in  this 
research  and  velocity  model  (I.N.2)  with  those  of  Anderson  and 
his  colleagues  at  the  California  Institute  cf Technology  (CIT  HOB). 
The  solid  line  curve  results  directly  from  observations  of  the 
travel  time  of  the  S  phase  for  a  surface  focus.  From  this  the 
I.N.2  velocity  model  is  deduced.  The  CIT  11GB  velocity  model, 
on  the  other  hand,  is  calculated  from  the  observed  dispersion  of 
surface  waves.  Then  from  this  velocity  model,  shown  in  the  figure 
by  dashed  lines,  the  dashed-line  travel-time  curve  Is  computed 
(Julian  and  Anderson,  1968).  The  overall  agreement  between  the 
two  travel-time  curves  is  considered  to  be  remarkably  good  except 
for  the  branches  AB  and  BC,  where  the  large  discrepancies  are 
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caused  by  a  different  depth  to  the  low- velocity  channel. 

It  should  be  noted  that  our  travel-time  curve  of  S  from 
nuclear  explosion  data  is  the  first  of  its  kind.  The  techniques 
which  make  the  identification  of  S  possible,  namely  resolution 
into  SH  and  SV  components  and  the  use  of  particle  motion  dia¬ 
grams,  were  developed  in  earlier  VELA-Uniform  work  done  by  us. 

S  wave  amplitude  data  from  the  Nevada  explosions  GREELEY  and 
HALF  BEAK  were  used  for  three  purposes,  namely:  l)  to  distinguish 
between  a  discontinuity  in  velocity  and  in  velocity  gradient  at 
depths  of  about  400  and  750  km  in  the  mantle,  2)  to  estimate  the 
magnitude  of  the  nuclear  explosions  from  the  S  wave  data,  and 
3)  to  study  the  pattern  of  radiation  of  seismic  energy  from  the 
explosions . 

Relatively  large  amplitudes  were  observed  near  point  E  of 
the  travel-time  curve,  but  there  was  no  increase  in  amplitude 
near  D  (see  Fig.  5).  This  suggests  that  the  velocity  increase 
beginning  at  750  km  depth  is  continuous,  but  that  the  velocity 
gradient  is  discontinuous  (Bullen,  196l).  There  is  no  indica¬ 
tion  of  an  amplitude  increase  near  points  B  and  C,  which  indi¬ 
cates  that  the  velocity  increase  at  400  km  depth  is  discontinuous. 

Differences  in  the  magnitudes  determined  from  P  wave  ampli¬ 
tudes  and  from  Rayleigh  wave  amplitudes  are  presently  considered 
to  be  an  important  discriminant  between  earthquakes  and  explosions. 
Additional  magnitude  estimates,  based  on  S  wave  amplitudes,  may 
be  of  help  in  such  discrimination.  From  our  amplitude  measure¬ 
ments  of  S  we  determined  the  magnitude  of  HALF  BEAK  to  be  4.9, 
with  a  range  of  values  of  4.5  to  5.2  and  a  standard  deviation 
of  0.2.  For  GREELEY  the  average  S  wave  magnitude  was  5.3,  with 
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a  range  of  4.8  to  3.8  and  a  standard  deviation  of  0.2.  Magnitude 
determinations  for  these  are  summarized  in  Table  1. 

Table  1 


Event 

P  wave 

Magnitude 

S  wave 

Surface  wave 

HALF  BEAK 

6.02  +  0.60 

(Teledyne) 

4.9  +  0.2 
(St.  Louis 

u.) 

6.2  (Uppsala) 

GREELEY 

6.29  +  0.45 
(Teledyne ) 

5.3  +  0.2 
(St.  Louis 

U.) 

5.6 

(Eskdalemuir) 

6.4  (Uppsala) 

6.1  (Eskdalemuir) 

Note  that  for  GREELEY  the  S  wave  magnitude  is  0.3  units  lower 
than  the  surface  wave  magnitude,  and  that  for  both  HALF  BEAK  and 
GREELEY  the  S  wave  magnitude  is  a  whole  unit  smaller  than  the  P 
wave  magnitude.  For  earthquakes,  on  the  other  hand,  the  P  and  S 
wave  magnitudes  are  the  same,  after  the  effects  of  the  radiation 
pattern  have  been  considered.  Thus  relatively  small  amplitude  S 
waves,  compared  to  P,  are  suggestive  of  an  explosion.  This  obser¬ 
vation  in  itself  is  of  not  much  importance  in  the  nuclear  explo¬ 
sion  detection  problem  because  only  very  large  explosions  generate 
S  waves.  But  the  converse  may  be  of  value,  i.e.,  any  event  for 
which  the  S  wave  magnitude  is  as  large  or  larger  than  the  P  wave 
magnitude  must  be  an  earthquake. 

The  radiation  patterns  for  both  HALF  BEAK  and  GREELEY  were 
similar  for  the  SH  arrivals  of  branches  AB  and  EF  and  the  hori¬ 
zontal  SV  arrivals  of  branch  EF,  which  taken  together  constitute 
most  of  the  S  arrivals.  The  SV  arrivals  were  all  positive 
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(towards  the  epicenter)  whereas  the  SH  arrivals  were  counter¬ 
clockwise  as  viewed  from  the  epicenter  for  stations  to  the  north¬ 
west,  north,  and  north-northeast  and  clockwise  for  stations  to  the 
east-northeast,  east,  and  east-southeast.  The  radiation  pattern 
of  SH,  considered  together  with  the  observation  that  the  first 
P  motion  appeared  as  a  compression  at  all  distances  and  azimuths. 
Indicates  that  the  focal  mechanism  can  be  represented  by  a  com¬ 
posite  of  a  uniform  compression  and  a  double  couple.  Toksoz  and 
Clermont  (1967),  using  surface  wave  and  P  wave  amplitudes,  found 
a  similar  kind  of  mechanism  for  BILBY,  If  one  assumes  that  the 
forces  of  the  double  couple  lie  in  a  horizontal  plane,  the  direc¬ 
tion  of  the  predicted  SH  nodal  line  for  BILBY  based  on  their  Ray¬ 
leigh  and  P  wave  study  agrees  with  that  observed  for  HALF  BEAK 
and  GREELEY.  This  leads  to  the  conclusion  that  the  radiation 
pattern  for  all  these  events  is  similar,  even  though  the  dis¬ 
tance  separating  the  shot  points  Is  as  much  as  fifty  km  and 
BILBY  was  in  tuff,  HALF  BEAK  in  rhyolite,  and  GREELEY  in  zeo- 
litized  tuff. 

dT/d  of  S  waves  across  LASA. 

Develecorder  records  of  the  long  period  seismographs  of  the 
D  and  F  rings  at  LASA  were  used  to  obtain  dT/d  values  for  the 
S  waves  across  the  array.  The  records  of  twelve  earthquakes, 
whose  hypocentral  coordinates  are  listed  in  Table  2,  were  digi¬ 
tized  and  then  resolved  into  SH  and  horizontal  SV  components  of 
motion.  Arrival  times  were  read  of  the  onset  and  the  first  peak 
or  trough  of  the  SH  motion.  The  waveforms  at  the  various  sta¬ 
tions  generally  were  consistent,  so  that  both  peaks  and  onsets 
gave  the  same  value  of  dT/d...  Particle  motion  diagrams  were 
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used  for  those  cases  in  which  the  onset  times  could  not  reliably 
be  determined  from  the  SH  component  directly. 


Table  2 


Origin 


Date 

time 

Lat 

• 

Long 

• 

Depth 

Magn 

1 

May 

67 

07-09-00.5 

39.7 

N 

21.3 

E 

15  km 

5 

.6 

9 

May 

67 

12-36-36.8 

56.6 

N 

152.6 

W 

33 

5 

.0 

11 

May 

67 

15-05-16.8 

20.3 

S 

68.5 

W 

67 

6 

.1 

13 

May 

67 

05-18-55.4 

56.5 

N 

152.6 

w 

33 

5 

.3 

16 

May 

67 

12-58-09.5 

13.5 

N 

90.6 

w 

95 

33 

4 

.8 

21 

May 

67 

07-18-12.8 

27.9 

N 

111.3 

w 

4 

.7 

27 

May 

67 

17-22-58.7 

51.9 

N 

176.1 

E 

34 

5 

.6 

19 

Jun 

67 

17-07-45.4 

52.7 

N 

166.9 

W 

33 

5 

.7 

15 

Nov 

67 

21-31-51.5 

28.7 

S 

71.2 

w 

15 

6 

.2 

17 

Nov 

67 

04-58-56.8 

28.5 

N 

43.8 

w 

33 

5 

.2 

23 

Nov 

67 

13-42-01.6 

80.2 

N 

1.0 

w 

10 

5 

.8 

1 

Dec 

67 

13-57-02.4 

49.5 

N 

154.4 

E 

136 

5 

.9 

The  values  of  dT/d;  ,  together  with  distance  to  the  center 
of  the  array,  are  presented  in  Table  3.  Also  included  are  the 
values  of  dT  and  dA  between  pairs  of  stations,  and  the  value  of 
dT/dA  from  the  Jeffreys-Bullen  tables  for  the  appropriate  dis¬ 
tance  and  focal  depth.  If  d  \  is  of  the  order  of  100  km,  an 
error  of  1  sec  in  dT  will  produce  an  error  of  about  0.01  sec/km 
in  dT/d\  .  The  latter  figure  corresponds  to  about  a  10#  error, 
because  for  teleseismic  distances  dT/d  .  varies  from  about  0.15 
to  0.07  sec/km.  This  example  illustrates  the  Importance  of  an 
accurate  determination  of  the  quantity  dT. 

Accurate  determinations  of  dT  require  an  adequate  sampling 
rate.  In  digitizing  the  Develocorder  records  we  sampled  approx¬ 
imately  four  points  per  second,  which  is  a  minimum  acceptable 
number  for  such  a  study.  One  point  per  second,  for  example,  is 
certainly  inadequate.  Because  the  waveforms  of  S  were  consistent 
across  the  array,  we  were  able  to  estimate  dT  to  the  nearest 
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0.1  sec  by  using  interpolation.  In  general,  the  most  reliable 
values  of  dT/d  _  are  those  for  which  d£  is  the  largest.  Only 
those  station  pairs  for  which  d  is  greater  than  100  km  are 
listed  in  Table  3» 


Table 


Station 

dT 

Date 

<°> 

Pair 

(sec ) 

1  May  67 

82.0° 

F1-F3 

18.1 

F3-F4 

10.7 

9  May  67 

29.9 

D2-F4 

16.8 

D4-F2 

18.1 

F2-F4 

27.1 

11  May  67 

74.8 

D2-F4 

12.9 

D4-F2 

14.2 

F2-F4 

21.5 

13  May  67 

29.9 

D2-F4 

17.2 

D4-F2 

18.1 

F2-F4 

27.5 

16  May  67 

35.6 

D2-F4 

17.6 

D4-F2 

16.8 

F2-F4 

26.6 

21  May  67 

19.1 

D1-F3 

21.9 

D2-F1 

21.5 

D3-F1 

21.1 

D4-F3 

22.3 

F1-F2 

31.4 

F1-F3 

35.7 

F2-F4 

27.1 

F3-F4 

31.4 

27  May  67 

48.7 

D2-F4 

14.6 

D4-F2 

16.3 

F2-F4 

24.1 

19  Jun  67 

38.7 

D2-F4 

16.3 

D4-F2 

15.9 

F2-F4 

25.8 

15  Nov  67 

81.4 

D2-F4 

14.9 

17  Nov  67 

51.3 

D3-F1 

11.1 

F1-F3 

18.3 

F1-F4 

17.1 

F2-F3 

16.2 

d  a 

dT/d 

dT/d  'avg 

dT/dj.B 

(km) 

(km/sec) 

(km/sec ) 

( km/sec ) 

211.89 

0.085 

132.73 

0.081 

0.083 

0.092 

116.97 

0.143 

124.17 

0.146 

189.22 

0.143 

0.144 

0.142 

123.46 

0.104 

133.40 

0.106 

200.72 

0.107 

0.106 

0.100 

116.84 

0.147 

123.98 

0.145 

188.98 

0.145 

0.146 

0.142 

122.25 

0.144 

135.53 

0.126 

199.44 

0.133 

0.138 

0.134 

114.91 

0.190 

109.54 

0.196 

113.21 

0.186 

116.88 

0.191 

162.45 

0.193 

185.95 

0.191 

135.39 

0.199 

158.89 

0.197 

0.193 

0.201 

116.87 

0.125 

123.86 

0.131 

188.90 

0.128 

0.128 

0.125 

113.38 

0.144 

119.66 

0.132 

183.09 

0.140 

0.139 

0.136 

123.16 

0.121 

0.121 

0.093 

101.30 

0.109 

150.98 

0.121 

132.49 

0.128 

123.32 

0.131 

0.122 


0.124 
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,6 . 

Station 

dT 

d_ 

dT/d  l. 

dT/d.:,  dT/d 

Date 

(  ) 

pair 

(sec) 

(km) 

( km/sec ) 

(km/sec)  (km/sec 

23  Nov  67 

46.9 

D3-F1 

F1-F2 

F2-F4 

14.0 

20.4 

17.5 

111.39 

163.44 

137.45 

0.125 

0.124 

0.127 

0.125  0.128 

1  Dec  67 

61.3 

D2-F4 

D4-F2 

F2-F4 

13.2 

14.1 

22.6 

120.74 

128.93 

195.58 

0.109 

0.109 

0.115 

0.111  0.113 

Comparison  of  the  observed  and  calculated  values  of  dT/d . 
shows  good  agreement  out  to  60°,  with  the  difference  between 
values  always  less  than  5 %.  There  are  an  insufficient  number 
of  data  to  look  for  any  azimuthal  dependence  such  as  has  been 
found  for  the  dT/d /  of  P  waves  across  LASA.  However,  we  plan 
to  continue  this  work  as  additional  usable  data  become  available 
and  hope  to  be  able  to  acquire  enough  data  to  investigate  the 
effect  of  azimuth. 

Spectral  Properties  of  Ground  Motion. 

The  theory  that  predicts  the  character  of  the  spectra  of 
the  components  of  the  S  wave  ground  motion,  which  depends  on 
earth  structure,  angle  of  incidence,  and  wave  frequency,  has 
been  developed  by  Haskell  (i960,  1962)  and  Nuttli  (1964).  The 
predicted  amplitude  vs  frequency  relation  for  SV  is  especially 
interesting  because  it  indicates  that  there  are  large  changes 
in  amplitude  with  small  changes  in  frequency  and/or  angle  of  in¬ 
cidence  for  particular  values  of  the  angle  of  incidence  and 
frequency.  However,  none  of  these  features  have  been  verified 
by  observation.  The  purpose  of  this  research  is  to  obtain  S 
wave  spectra  from  the  observational  data  and  to  compare  some  of 
their  properties  with  those  expected  from  theory. 

Seismograms  of  five  earthquakes  have  been  analyzed.  Table  4 
lists  their  dates  and  hypocentral  coordinates.  The  long  period 
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Table  4 


Earth¬ 

quake 

Number 

Date 

Origin 

Time 

Ld  t  1 

Long. 

Depth 

Magn. 

1 

13  Apr  63 

02-20-37.9 

6.26°S 

76. 66°W 

125  km 

6  3/4-7 

2 

15  Oct  63 

09-59-30.1 

67.20  N 

18.40  W 

33 

3 

12  Jan  64 

06-00-13.2 

53.20  N 

166.30  W 

33 

5§ 

4 

28  Jun  64 

17-27-59.8 

4.00  N 

32.40  W 

33 

5.3 

5 

18  Sep  64 

13-12-42.3 

39.80  N 

29.70  W 

20 

seismograms  used  are  from  LRSM  stations  at  eplcentral  distances 
o  o 

of  50  to  70  .  Table  5  shows  the  stations  used  for  each  event, 
their  eplcentral  distances,  and  back  azimuths. 


Table  5 

Earthquake 

Azimuth  (station 

Number 

Station 

Distance 

to  epicenter) 

1 

DHNY 

48.23° 

182.15° 

DRCO 

52.33 

139.29 

FMUT 

55.97 

135.56 

LCNM 

47.90 

137.78 

NOWS 

52.76 

165.39 

PMWY 

54.08 

143.66 

SEMN 

52.97 

157.15 

TFCL 

58.01 

126.25 

WNSD 

53.70 

150.28 

2 

AZTX 

55.44 

28.05 

BLWV 

45.36 

29.19 

RTNM 

55.22 

28.25 

SKTX 

55.14 

28.06 

3 

DRCO 

56.62 

68.77 

FRMA 

52.99 

75.25 

GIMA 

51.24 

77.40 

JUTX 

54.24 

66.34 

4 

BRPA 

55.04 

118.11 

BLVJV 

55.83 

114.66 

DRCO 

75.99 

95.76 

LSNH 

53.04 

127.39 

5 

DRCO 

59.15 

61.40 

JELA 

50.28 

62.49 

RTNM 

57.09 

62.25 

RYND 

50.62 

71.22 

WNSD 

51.45 

68.22 

In  order 

to  obtain  the 

spectral 

components  of  S  we  em 

ployed  the  following  procedure.  The  seismograms  were  digitized 


at  an  0.5  sec  interval,  and  the  horizontal  components  of  motion 
resolved  into  components  along  (UR)  and  transverse  to  (UT)  the 
great  circle  path  at  the  station.  The  average  length  of  seismo¬ 
gram  digitized  corresponded  to  about  forty  seconds  of  time.  Be¬ 
cause  S  is  followed  closely  by  other  phases  such  as  sS  and  PS, 
the  data  were  passed  through  a  Hamming  lag  window  to  suppress 
the  effect  of  these  later  phases.  After  reductions  for  the  fre¬ 
quency  response  of  the  seismograph  system  were  made,  the  data 
were  Fourier  analyzed  using  a  technique  known  as  Filon's  method. 

The  results  to  date  of  the  analysis  can  best  be  presented 
in  the  form  of  some  typical  examples.  Theoretical  curves  of  the 
spectral  ratio  FR/FZ  were  prepared  for  the  Haskell  and  McEvilly 
crustal  models  (see  Table  6  for  properties)  for  angles  of  inci¬ 
dence  between  20°  and  36°  in  2°  steps.  Curves  of  these  spectral 
ratios  were  also  prepared  for  a  modified  Haskell  model,  in  which 
the  crustal  velocities  were  changed  to  =*  =  6.6  km/sec  and 
fi  »  3.86  km/sec. 

Table  6 

Comp,  wave  Shear  wave 


Model 

Layer  thick¬ 
ness  (km) 

velocity 

(km/sec) 

velocity 
(km/sec ) 

Density 

(gm/cm3) 

Haskell 

37.0 

6.285 

3.635 

2.869 

infinite 

7.960 

4.600 

3.370 

McEvilly 

1.0 

4.40 

2.50 

2.70 

10.0 

6.10 

3.50 

2.70 

9.0 

6.30 

3.66 

2.80 

18.0 

6.60 

3.86 

2.90 

infinite 

8.15 

4.70 

3.30 

In  discussing  the  plots  of  the  spectral  ratios  we  shall 
use  the  terms  "character"  and  "shift"  as  defined  by  Hannon  (1964). 
"Character"  refers  to  the  general  shape  of  adjacent  peaks  and 
troughs  of  the  spectral  ratio  curves.  One  curve  is  "shifted" 
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with  respect  to  the  other  when  segments  of  the  curves  having  sim¬ 
ilar  character  are  displaced  with  respect  to  the  frequency  co¬ 
ordinate  . 

Figure  6  presents  the  ratio  FR/FZ  of  the  amplitude  spectra 
for  earthquake  number  1  at  the  station  FMUT,  as  well  as  a  theo¬ 
retical  curve  for  McEvilly's  earth  model  for  an  angle  of  inci¬ 
dence  of  3^°.  The  left  part  of  both  curves  (frequencies  less  than 
0.18  cps)  shows  the  same  character.  The  difference  in  amplitudes 
of  the  curves  can  be  reduced  by  decreasing  the  velocity  contrast 
between  the  half  space  and  the  lowest  crustal  layer  in  the  theo¬ 
retical  model.  Better  agreement  between  the  theoretical  and 
observed  curves  also  might  be  obtained  by  a  slight  increase  in 
total  crustal  thickness  of  the  theoretical  model. 

Observed  and  theoretical  spectral  ratio  curves  of  FR/FZ 
at  TFCL,  also  for  earthquake  number  1,  are  given  in  Fig.  7.  The 
angle  of  incidence,  34°,  is  the  same  as  for  Fig.  6,  which  is  ex¬ 
pected  because  both  stations  are  at  about  the  same  epicentral 
distance.  The  modified  Haskell's  model  gives  a  theoretical 
curve  which  ha3  the  same  character  as  the  observed  one.  The 
different  levels  of  amplitude  indicate  that  further  adjustments 
are  still  required  in  the  velocities  of  the  theoretical  model. 

A  shift  of  the  theoretical  curve  to  the  left  is  also  needed. 

Figure  8  contains  spectral  ratios  of  the  station  BRPA  for 

earthquake  number  4.  The  best  fitting  theoretical  curve  is  for 

o 

an  angle  of  incidence  of  32  ,  less  than  that  of  Figs.  6  and  7, 
which  can  be  explained  by  the  fact  that  the  focal  depth  of  earth¬ 
quake  number  1  is  the  greater.  A  change  in  amplitude,  which  can 
be  produced  by  changing  the  velocities  of  the  model,  is  required. 
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The  observed  spectral  ratio  curve  of  LSNH  for  earthquake  number 
4  (Pig.  9)  has  the  position  of  the  peaks  and  troughs  Interchanged 
as  compared  to  those  of  Pig.  8.  This  Is  surprising,  because  both 
stations  are  approximately  at  the  same  distance  and  might  be  ex¬ 
pected  to  have  the  same  crustal  structure. 

Figure  10  Is  unusual  In  the  sense  that  the  best  fit  for 
RTNM  for  earthquake  number  5  Is  found  to  be  McEvilly's  model  with 
an  angle  of  Incidence  of  only  20°.  The  agreement  In  character 
of  the  observed  and  theoretical  curves  Is  fairly  good,  although 
the  theoretical  curve  should  be  shifted  slightly  to  the  left. 

In  the  previous  discussion  the  best  fitting  theoretical 
curve  was  taken  to  be  one  of  a  set  of  curves  for  the  McEvllly 
or  Haskell  models  at  angles  of  Incidence  of  20,  22,  24,  26,  28, 
30,  32,  34,  35  and  36  degrees.  It  la  apparent  that  none  of  these 
theoretical  curves  gives  an  Ideal  match  to  any  of  the  observed 
curves.  Better  agreement  can  be  had  by  changing  the  parameters 
of  the  earth  model  used  to  compute  the  theoretical  curves.  Fur¬ 
ther  effort  will  be  devoted  to  finding  a  systematic  way  of  chang¬ 
ing  the  parameters  so  as  to  get  a  better  fit,  such  as  by  using 
partial  derivatives  of  the  spectral  ratio  with  respect  to  the 
various  parameters . 
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4.  Earthquake  Energy  Determination 
The  present  method  of  energy  determination  from  the  magni¬ 
tude-energy  relations  of  Gutenberg  and  Richter  or  others  is  not 
very  accurate.  With  the  advent  of  the  modern  high-speed  elec¬ 
tronic  computers,  with  the  better  understanding  of  focal  mechan¬ 
ism  of  earthquakes  and  Q  of  the  earth,  with  the  knowledge  of 
well-determined  crustal  structure  and  the  existence  of  highly 
sensitive  standard  seismographs  with  we 11 -determined  instru¬ 
mental  response  characteristics,  it  is  now  possible  for  us  to 
take  into  account  the  various  factors  affecting  the  ground  mo¬ 
tion  and  to  come  up  with  more  accurate  energy  estimates  than 
it  was  possible  in  Gutenberg's  time. 

The  seismogram  can  be  visualized  as  the  output  of  a  series 
of  filters  into  which  an  original  seismic  impulse  has  been  fed. 
Each  filter  corresponds  to  a  portion  of  the  propagating  medium 
and  is  characterized  by  a  complex  transfer  function  H(^),  which 
operates  on  the  input  spectrum  and  modifies  it  as  it  passes 
through.  The  final  output  of  this  filter  series  is  P  (<•'),  which 
is  measurable  at  the  surface.  These  complex  functions  may  be 
written  in  terms  of  real  amplitude  and  phase  spectrums: 


(1) 

(2) 

(3) 
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S(^)  -  A0(ca)  e1  & 

F(m)  =  A(^)el7i(^ 
HDM  m  AD(^)eifD(a') 


and  similarly  for  HCR(<o)  and  H  INS(^)*  where  H^),  hcr(^), 
HiNS(m)  are  the  complex  transfer  functions  of  the  mantle,  crust 
and  the  seismograph.  It  then  follows  from  the  definition  of 
the  transfer  function  that 


Aq(«)  ~  i  AINS^'')  ACR^’)  AD^^j  (*0 

^0(>-')  ”  +  *CR<'a>  ^DM.l  (5) 


where  A  (o )  and  f  ('■>)  are  known  from  the  observed  body  wave  data. 

A  program  has  been  written  for  obtaining  discrete  Fourier 
transforms  from  the  digitized  data  of  the  pulse  recorded  on  the 
seismogram.  An  existing  program  has  been  adapted  for  obtaining 
the  steady  state  harmonic  response  of  the  seismograph.  There  is 
provision  for  introducing  magnification  values  at  five  frequen¬ 
cies  in  the  input  data  and  the  response  for  any  intermediate 
frequency  is  obtained  by  linear  interpolation.  A  program  based 

on  the  Haskell-Thomson  matrix  formulation  for  obtaining  the  crustal 
response  for  incident  P  and  S  waves,  as  a  function  of  frequency 
and  angle  of  Incidence,  was  tested. 

As  in  Teng  and  Ben-Menahem  ( 1965 ) ,  the  amplitude  response 
for  the  effect  of  geometrical  spreading  and  absorption  in  the 
mantle  may  be  expressed  as 

Ap(u)  »  G  x  AT(^)  (6) 


°in  ^  dihT/:> 

cos  iQ  sin  ©  d©  J 


(7) 


where 
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is  a  factor  due  to  geometrical  spreading  and  AT(ci)  is  the  atten¬ 
uation  function.  Assuming  a  Bullen  A  velocity  model  and  a  Q(H) 
model  calculated  by  Anderson  and  Archambeau  (1964),  Teng  and  Ben- 
Menahem  obtain  the  expression  for  the  attenuation  function  as 
follows: 


AT(.  )  -  exp[-jray 

„e4-,^131n(l^0>lxi°-3] 


(8) 


where  the  subscripts  o  and  h  refer  to  the  base  of  the  crust  and 
the  source  level,  respectively.  A  program  has  been  prepared  to 
correct  for  the  effect  of  geometrical  spreading  and  absorption 
in  the  mantle  and  reduce  it  to  the  focal  sphere. 

Let  (x,y,z)  be  the  coordinates  of  a  point  in  a  right-handed 
coordinate  system  X,  Y,  Z  related  to  the  orientation  of  the 
source  mechanism,  let  (x,y,”z)  be  the  coordinates  of  the  same 
point  in  a  geographic  coordinate  system  (N,  E,  and  down)  re¬ 
ferred  to  the  origin  of  the  focal  sphere,  and  let  (oV^x'  'V 
and  («  y,/&y,  Y~y)  be  the  direction  cosines  of  the  X  and  Y  axes 
with  respect  to  x,  y,  z.  Then  the  displacement  field  due  to  a 
double  couple  source  with  forces  oriented  along  the  X  and  Y 
axes  is  given  by  (see  e.g.,  Bessonova  et  al,  i960) 

(9) 


ua  .  2  *  1L,  k'  (t-fi) 
a  I?  r3  a' 

k'(t-H) 


H 


P 

ub 


_  1 _ _  C 

^rrp  Pr 

1  D 
^TTp  pp 


*’  <t-§) 


(10) 

(n) 


H  P 

where  u.  is  the  radial  motion  of  the  P  wave,  ub  and  ub  are  the 


SH  and  SV  components,  respectively,  a  and  b  are  the  velocities 
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of  P  and  S  waves,  respectively,  R  is  the  radial  distance  to  the 
point  of  observation  and 


(  <X  y  +  ^  y  x)  y  -(/$£_ 

y  +/*y  x)  X 

(12) 

sin  1 

o 

=  2xy  cos21q  -  z  (/*x  y  + 

fy  X) 

(13) 

sin  10  cos  10 


Based  on  these  formulae  a  program  has  been  prepared  for  the 
P  and  S  wave  radiation  pattern  on  the  focal  sphere  as  a  function 
of  the  angle  of  Incidence  at  the  focus  and  azimuth,  due  to  a 
double  couple  source  of  arbitrary  orientation. 

The  next  step  is  to  develop  a  formula  for  the  energy  deter¬ 
mination  from  the  observed  spectrum  of  body  wave  pulses  and  to 
combine  the  above  programs.  At  a  later  stage  the  effect  of  the 
finiteness  of  the  source  and  rupture  velocity  is  planned  to  be 
incorporated  in  the  energy  formula.  At  present  attention  is  con¬ 
centrated  on  intermediate  and  deep  focus  earthquakes  which  give  P 
and  S  pulses  isolated  from  pP  and  sS,  respectively.  Later  shallow 
earthquakes  may  be  included.  For  these  the  effect  of  interference 
with  the  main  pulse  of  crustal  reverberations  must  be  removed 
prior  to  analysis. 
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5.  The  Effect  of  Focal  Depth  on  the  Spectrum  of  Body  Waves. 

As  an  approach  to  the  problem  of  determining  the  focal 
depth  of  a  seismic  event  from  the  spectrum  of  body  waves,  the 
following  problem  was  considered:  In  a  layered  spherical 
earth  a  system  of  forces  acts  at  a  particular  level.  The  prob¬ 
lem  is  to  determine  how  the  ratio  of  the  frequency  spectrum  of 
body  waves  recorded  at  a  point  on  the  surface  of  the  sphere  to 
the  frequency  spectrum  of  the  system  of  forces  at  the  source 
level  is  influenced  by  the  depth  of  the  layer  at  which  the 
force  system  acts.  This  problem  was  investigated  for  a  layered 
sphere  (having  arbitrary  number  of  layers  and  a  liquid  core) 
with  arbitrary  orientation  of  the  system  of  forces,  acting  at 
one  particular  level.  An  analytical  solution  using  Haskell- 
Thomson  type  of  complex  matrices  was  obtained.  However,  the 
numerical  calculation  even  for  the  simplest  kind  of  force  dis¬ 
tribution  appeared  very  difficult,  and  the  spherical  earth 
model  was  put  a3ide.  The  problem  was  reformulated  for  a 
layered  flat  earth  model. 

The  system  consists  of  (n-1)  layers  over  a  half-space. 

The  free  surface  is  labelled  "0",  and  the  ith  interface  is 
the  bottom  of  the  ith  layer.  The  half-space  is  layer  n.  On 
the  3  interface  a  system  of  time-varying  forces  acts.  The 
problem  is  to  derive  the  amplitude  of  plane  waves  which  propa¬ 
gate  along  a  specified  direction  to  a  point  in  the  interior  of 
the  system.  This  direction  has  a  specified  azimuth  with  respect 
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to  the  epicenter  and  a  specified  inclination  with  respect  to 
the  vertical  at  the  focus. 

Inside  the  half-space,  at  distances  far  away  from  the 
source  layer,  the  wave  propatation  along  the  specified  direc¬ 
tion  can  be  considered  as  equivalent  to  wave  propagation  along 
a  ray  path  which  is  characterized  by  a  particular  azimuth  and 
a  particular  ray  parameter.  Thus,  the  amplitude  of  plane  waves 
of  various  frequencies  which  propagate  in  the  half-space  will 
give  considerable  information  regarding  the  spectrum  of  body 
waves  which  one  expects  to  record  at  a  particular  epicentral 
distance  and  at  a  particular  azimuth  for  a  spherical  earth. 

In  particular,  it  is  possible  to  investigate  how  the  amplitudes 
of  plane  waves  of  various  frequencies  are  influenced  by  the 
depth  of  the  layer  on  which  the  system  of  forces  act.  The 
formulation  of  the  problem  follows. 

Mathematical  formulation. 

Displacement  potentials  are  defined  as 

*$■  O'  v'  w)  =  9ra<*  9  **  Ca^  CUY'  *  cuY^  (.0>0'  )  ^ 


Equatiore  of  motion  are  satisfied  if 
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^  d-t 


The  displacements  and  stress  are  given  by 
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where  f1,  fg,  f^,  are  functions  of  <p  and  y  and  and  fg  are 
functions  of  A  alone: 
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The  potential  A  does  not  contribute  to  the  ver¬ 
tical  components  of  displacement  and  normal  stress  and  thus 
represents  SH  waves.  The  potential  <p  represents  P  waves  and 
y  represents  SV  wave  : .  P  and  SV  waves  are  considered  together 
and  SH  waves  separately. 

P  and  SV  Waves. 

We  consider  plane  wave  solutions  of  equation  (2)  in  the 
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The  displacements  and  the  stresses  will  be  continuous  across 
the  various  interfaces  if  the  functions  f1#  f2,  fi±  are  con 
tinuous  across  these  interfaces.  Using  eq.  (4),  the  functions 
fl  ...  f^  can  be  written  as: 
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Let  us  consider  a  surface  Z  »  constant.  On  this  surface 

fi  (x,  y,  z)  can  be  written  as  a  two-dimensional 

Z  =  const. 

Fourier  expansion  as  follows: 

[  Fi  (*.**)  e‘  ^K'1C  dK,  ^K, 
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Functions  F^  (^K^)  are  functions  which  are  defined  on  the  sur- 
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faces  Z  =  constant.  Prom  eqs.  (7)  and  (6),  one  may  write  the 
matrix  relation 


The  elements  of  the  matrix  R  are: 
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Now  consider  the 
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layer.  It's  upper  surface  is  (m-l)th 
Placing  the  origin  of  the  coordinate  system  at 
interface,  using  eq.  (9),  one  can  write 
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where  the  elements  of  the  matrix  B  are: 
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Let  the  thickness  of  the  mth  layer  be  Then  using 


one  can  write 


where  the  elements  of  the  matrix  A  are 
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From  eqs.  (10)  ami  (IS) 
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<“> 

eq.  (9). 
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Prom  eq.  (10)(*) 
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Using  eq.  (14),  we  have 
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Let  the  4  x  4  matrix  Bn”1  An-1 
be  denoted  as  the  matrix  P . 

Using  eq.  (13)*  one  obtains 
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*  Since  there  are  no  source  in  the  half-space,  *^2  “ 

(  4)2n  #  y2n  represent  amplitude  of  plane  P  and  SV  waves  respec¬ 
tively  which  travel  upward  in  the  half-space) , 


where  F^,  F2°,  etc.,  refers  F1(K1K2)  on  the  0th  interface, 
i.e.  on  the  free  surface. 

We  consider  that  at  the  source  level  the  displacements 
are  continuous  but  the  stresses  are  not.  Let  Ap£2(x,y) 
the  discontinuity  in  the  normal  stress  and. let  the  discontin¬ 
uity  in  tangential  stress  PZx(x,y)  be -^-(h(x,y))  and  that  in 
the  tangential  stress  Pgy(x*y)  be  l y )).  Let  ^(K^Kg) 
and  be  the  two-dimensional  Fourier  transform  of  h(;t,y) 

and  Apzz(x,y)  respectively.  Then 


Therefore,  using  eq.  (18),  one  obtains 


Let  the  matrix  product  ....  B8A8  be  denoted  as  the 

matrix  Q. 

Using  eqs.  (17)  and  (19)>  one  obtains 


Note:  Sth  interface  is  S+  on  the  downward  side  and  S"  on  the 
upward  side. 
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Solving  eq.  (20)  for  the  four  unknowns  i1!94*  ?2S+>  P3S+> 
one  can  obtain  from  eq.  (17) 
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4*2°  (K^Kg)  and  H'1n(K1Kg)  are  the  amplitude  of  P  and  SV  waves 
respectively  in  the  half-space. 

SH  Waves. 

Displacements  and  stresses  are  given  by 

?-  is  <.*•?.’> 

-JU  *  <-»  •  •V'J  (’-'i 

v«i  *  ^-i-v 

Functions  f-  (x,y, 2),  fg  (x,y,z)  are  given  following  eq.  (3). 
Solution  of  wave  eq.  (2)  is 
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U8ing  eq.  (22),  functions  f^(x,y,z),  fg(x,y,z)  can  be  written 
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On  the  surface  Z  ■  constant,  fc  (x,y,z)\  .  and 

D  1Z  *  cons u . 

tc  (x,y,z)i  can  be  written  as  two-dimensional  Fourier 

u  i Z  w  const. 

Integral  expansion  as  follows: 
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Prom  cqs.  (23)  and  (24)  one  can  write 
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Let  the  2x2  matrix  on  the  left-hand  side  of  eq.  (25)  be 
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denoted  as  L.  All  the  elements  of  the  matrix  L  are  to  be 
evaluated  at  Z  *  constant. 

Now  consider  the  mth  layer  again.  Placing  the  origin  of  the 
coordinate  at  (m-l)*'11  Interface,  using  eq.  (25),  one  obtains 
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where  the  elements  of  the  matrix  D  are: 
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Let  the  thickness  of  the  rath  layer  be  d  ,  Then,  using  eq.  (25) 
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where  the  elements  of  the  matrix  C  are 
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From  eqs.  (26)  and  (27),  one  can  write 
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Prom  eq.  (26),  one  obtains 
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Using  eq.  (29),  one  obtains 
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Let  the  matrix  product  Dn  Cn.]_  ^n-1  ' ’ ’  s+^  3 

denoted  as  the  matrix  X. 

Using  eq.  (28),  one  obtains 
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Two-dimensional  Fourier  transform  of  different  kinds 
stress  distributions 

H(K}K2)  functions  for  P  and  SV  waves. 

(1)  a  single  vertical  force  acting  at  she  origin  in  the  posi¬ 
tive  Z-direction: 
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(2)  a  single  vertical  couple  with  counterclockwise  moment 
about  an  axis  in  the  x-y  plane,  making  angle  with  the 
x-axls  -  the  moment  arm  is  . 
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(3a)  two  couples  with  opposing  senses  in  the  x-y  plane.  The 
moment  arm, It  ,  makes  an  angle  with  the  X-axis 
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taking  Fourier  transform  of  both  sides,  one  obtains 

where 

^■2>,  ^  are  the  Fourier  transform  of  \yK  C<-'0 
and  ^  respectively.  It  turnG  out  that 
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The  above  system  reduces  to  a  single  couple  if  either 
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.  In  that  case  in  the 
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expression  for  H(KjK2)  above  one  can  neglect  either  Mgx  or  MZy* 
In  case  »  -Mzx  -  M,  then 
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(3b)  horizontal  double  couple 
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taking  two-dimensional  Fourier  transform  of  both  sides,  one 
obtains 

- r («’-**')  »<*'**)  *  *'*'**y 

where  and  FZy  are  Fourier  transform  of  pZx(x,y)  and 
pZy(x,y)  respectively.  It  turn  out  that 
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HfK-^)  functions  for  SH  vraves. 

A  pair  of  horizontal  couples. 
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taking  Fourier  transform  of  both  sides,  one  obtains 
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Discussion 

A  system  of  forces  of  harmonic  time  variation  having 
arbitrary  orientation  acts  on  the  S  interface.  The  amplitude 
of  plane  P,  SV  and  SH  waves  which  propagate  along  a  direction 


fn • the  half-space  are  expressed  In  terms  of  the  two-dimensional 
spatial  Fourier  transform  HfK^Kg)  of  force  distribution  on  the 
source  level.  Two-dimensional  Fourier  transform  H(K1,H2)  is 
calculated  for  different  kinds  of  stress  distribution.  If  these 
functions  are  used  In  eq.  (20)  (in  the  case  of  P  and  SV 
waves)  and  In  eq.  (33)  (in  the  case  of  SH  waves),  one  would  ob¬ 
tain  by  using  eqs.  (20a)  and  (33a)  the  amplitude  of  P,  SV  and 
SH  waves  which  propagate  along  a  particular  direction  In  the 
half-space. 

At  the  present  time  computer  programs  are  being  written 
for  various  phases  of  computational  work.  The  object  of  com¬ 
putation  is  to  calculate  <$>in,  Yin,  A  1*1  for  different  w  for 
various  distributions  of  forces  at  the  Bource  level  and  various 
depths  of  the  source  interface. 
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6.  Model  Studies 
Analogue  Spectral  Techniques 

The  research  In  seismic  model  techniques  originated  during 
the  "period  of  this  report  was  devoted  primarily  to  the  develop¬ 
ment  of  rapid  analogue  techniques  for  extracting  spectral  Informa¬ 
tion  from  the  signals  generated  In  the  model  system.  In  previ¬ 
ous  model  Investigations  of  wave  generation  and  propagation, using 
primarily  explosive  sources,  but  also  pulsed  ceramic  transducers, 
analysis  was  done  In  the  time  domain  or  In  the  frequency  domain, 
using  digital  methods  to  extract  the  spectrums. 

The  acquisition  of  two  new  Items  of  laboratory  Instrumenta¬ 
tion  by  the  Department  has  made  possible  this  development  of  the 
capability  of  the  seismic  models  laboratory.  One  Item,  suitable 
for  on-line  examination  of  the  main  spectral  characteristics  of 
the  signal  produced  by  a  repetitively-pulsed  source.  Is  a  Nelson- 
Ross  Model  013  spectrum  analyzer.  The  other  Is  a  Kay  Electric 
Company  Model  7029A  Sound  Spectrograph. 

The  Nelson-Ross  analyzer  plugs  Into  a  Tektronix  535A  Oscil¬ 
loscope.  It  can  ba  used  to  display  the  entire  spectrum  of  a 
repetitive  signal  from  0  to  500  kHz,  or  any  desired  frequency  can 
be  selected  as  a  center  frequency  of  the  display  and  by  means  of 
a  dispersion  control  a  detailed  presentation  of  any  portion  of 
the  spectrum  can  be  obtained.  An  auxiliary  function  generator 
Is  used  as  a  control  on  the  calibration.  The  analyzer  requires 
a  repeated  signal,  and  so  far  has  only  been  applied  to  pulsed 
transducers  and  continuous  wave  inputs.  The  system  can  be  applied 
to  explosion  signals  by  recording  on  magnetic  tape  and  forming 
a  tape  loop. 


A  disadvantage  of  the  present  method  of  application  18  that 
the  entire  signal.  Including  all  body  and  surface  waves.  Is  pro¬ 
cessed.  It  Is  planned  to  Incorporate  a  gating  circuit  that  will 
pass  only  a  preselected  portion  of  the  signal,  e.g,,  the  first 
arriving  P-wave. 

The  first  application  of  the  Kay  Sound  Spectrograph  to  the 
analysis  of  seismic  transients  was  reported  by  Ewing,  et  al. 
(1959)  and  Landlsman,  et  al,  (1962).  In  their  work  It  was  neces¬ 
sary  to  speed  up  the  playback  of  the  data  In  order  to  put  the 
signal  Into  the  range  of  the  spectrograph.  In  model  seismology 
the  problem  Is  the  reverse:  the  energy  Is  mostly  In  the  range 
of  tens  of  kllo-Hertz,  so  that  a  slowed-down  replay  Is  required. 
This  has  been  achieved  by  recording  the  signals  on  a  Precision 
Instrument  Model  PS-207A  Instrumentation  Recorder  at  60  lps, 
direct  record,  and  playing  back  at  1  7/8  lps.  Frequencies  are 
divided  by  32,  putting  the  signal  Into  the  0-16  kHz  range  of 
the  spectrograph.  The  chief  limitation  is  the  100  kHz  high 
frequency  cut-off  of  the  tape  recorder. 

The  two  displays  available  from  the  spectrograph  that  are 
most  useful  for  this  work  are  Illustrated  In  Figures  11  and  12. 
Figure  11  Illustrates  the  Sonogram,  a  frequency-amplitude-time 
display.  The  horizontal  axis  Is  time,  the  vertical  axis  Is  fre¬ 
quency  increasing  upward,  and  the  amplitude  is  shown  by  relative 
darkness.  Frequencies  Indicated  next  to  the  calibration  marks 
have  been  converted  to  real  frequencies  and  the  time  scale  cor¬ 
responds  to  original  real  time.  The  arrival  of  the  transient 
from  the  explosion  In  Plexiglas  Is  clearly  seen,  followed  by 
low  frequency  reverberations.  Although  the  time  scale  is  very 
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compressed,  a  detailed  interpretation  is  possible  by  the  tech¬ 
nique  discussed  next. 

The  displays  in  Figure  12  are  called  sections.  A  section 
can  be  taken  at  any  preselected  point  in  the  signal.  Amplitude 
versus  frequency  (increasing  downward)  is  shown  at  the  times 
selected.  By  moving  through  the  signal  in  small  increments  of 
time,  a  complete  picture  of  the  growth  of  the  spectrum  as  time 
passes  can  be  assembled.  The  analysis  of  the  data  in  Figure  12 
is  discussed  below. 

A  variety  of  ests  were  made  to  determine  the  effects  on 
the  signal  of  processing  through  several  instruments.  The  output 
of  the  detectors  as  developed  on  an  oscilloscope  was  photographed 
and  used  as  the "true "record  for  comparison  purposes.  The  signal 
as  recorded  both  in  the  tape  recorder  and  in  the  tape  loop  of  the 
spectrograph  was  played  into  a  Brush  oscillograph.  It  was  clear 
that,  in  the  lower  frequencies  at  least,  little  is  lost  in  going 
from  the  recorder  to  the  spectrograph.  The  tape  recorder  itself 
gives  a  faithful  reproduction  when  played  back  at  the  slow  speed. 
This  test  was  of  limited  value  because  the  Brush  oscillograph 
itself  did  not  pass  signals  above  2.5  kHz,  so  that  the  fidelity 
at  the  higher  frequencies  could  not  be  tested. 

A  measure  of  the  reproducibility  may  be  seen  in  the  seismo¬ 
grams  in  Figure  13.  The  records  for  13a  and  b  are  from  two  differ¬ 
ent  shots  in  Plexiglas.  In  13c,  the  signal  that  produced  13b  has 
been  put  through  the  magnetic  tape  recorder  and  played  back  at 
1/32  of  the  signal  speed.  The  time  scale  on  the  oscilloscope 
has  been  adjusted  accordingly.  Unfortunately,  because  of  trigger¬ 
ing  problems,  the  onset,  the  most  interesting  part,  is  not  seen 


66 


In  13c,  but  the  rest  of  the  record  shows  that  the  signal  Is  faith¬ 
fully  reproduced  in  the  slow-down  procedure. 

A  thorough  calibration  of  the  spectrograph  was  conducted, 
using  continuous  sine  waves  and  square  wave  pulses.  Both  the 
frequency  calibration  and  the  time  calibration  along  the  record 
were  checked. 

As  an  Independent  test,  the  overall  frequency  response  of  the 
entire  system,  source  through  medium  and  receiver,  was  determined 
by  driving  the  source  with  a  sinusoidal  voltage  with  variable 
frequency.  Peaks  in  the  output  were  noted  and  compared  vlth 
those  determined  by  the  two  spectrum  analyzer*.  The  most  impor¬ 
tant  finding  of  this  series  of  experiments  was  that  most  of  the 
spectral  peaks  were  associated  with  sharp  resonances  in  the  driv¬ 
ing  and  receiving  transducers.  Various  mountings  and  transducer 
holders  were  tried,  but  some  strong  resonances  still  occurred 
within  the  part  of  the  spectrum  in  which  the  pulse  employed  for 
generating  transients  had  most  of  its  energy.  In  addition,  the 
spectrum  proved  sensitive  to  the  pressure  with  which  the  crystal 
was  held  against  the  medium.  For  those  reasons  explosive  sources 
and  capacitance  transducers  were  used  in  subsequent  studies .  This 
source  and  receiver  system  has  been  proven  through  many  years  of 
work  in  this  laboratory.  The  contrast  in  the  signals  produced  by 
crystal  and  capacitance  detectors  may  be  seen  in  Figures  14a  and 
14b,  respectively.  The  source  is  an  explosion  in  Plexiglas.  The 
capacitor  transducers  are  less  sensitive  than  the  ceramic  ones, 
so  that  the  voltage  from  the  function  generator  was  too  low  to 
produce  a  detectable  signal. 
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Preliminary  Experiments. 

As  a  test  of  the  methods  of  analysis,  a  series  of  experi¬ 
ments  on  the  frequency  dependence  of  transmission  through  a  layer 
over  a  half-space  was  conducted.  This  is  a  problem  receiving 
much  attention  in  seismology  today  in  an  effort  to  account  for 
the  effects  of  crustal  structure  on  the  recorded  surface  motion. 
When  the  effectiveness  of  the  analogue  techniques  and  the  exper¬ 
imental  methods  are  firmly  established,  related  problems  of 
greater  mathematical  difficulty,  such  as  transmission  through 
dipping  layers  will  be  attacked. 

In  the  experiments,  the  source  was  far  enough  from  the  sur¬ 
face  (five  wavelengths  or  more)  that  plane  wave  theory  was  appli¬ 
cable.  Only  normal  incidence  at  the  base  of  the  layer  was  con¬ 
sidered. 

In  the  first  series  of  experiments,  aluminum  was  used  as 
the  half-space,  with  Plexiglas  as  the  layer.  First,  the  spec- 
trums  were  observed  at  a  point  in  the  interior  of  the  half-space 
and  at  the  free  edge  (before  the  Plexiglas  layer  was  put  in  place) 
for  a  lead  titanate-zirconate  source  driven  by  both  a  repeated 
pulse  and  a  sine  wave  generator.  Then  the  layer  was  added  and 
the  observations  on  the  top  of  the  layer  made. 

Because  the  input  signal  in  aluminum  for  a  pulse  source  is 
rich  in  high  frequencies,  the  magnetic  tape  recorder  was  in¬ 
adequate.  Spectrums  obtained  with  the  Nelson-Ross  analyzer  and 
by  observing  the  response  to  continuous  sine  waves  are  shown  in 
Figures  15  and  16.  In  both  figures  ceramic  transducers  were 
used  as  both  source  and  receiver.  The  source  was  in  the  inter?. or 
of  an  aluminum  sheet.  The  receiver  was  at  the  surface.  The 
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results  show  that  the  analyzer  gives  a  satisfactory  picture  of  the 
spectrum.  The  peaks  in  these  figures  correspond  to  resonances  in 
the  transducers . 

It  was  known  from  previous  experience  that  impulsive  sources 
in  Plexiglas  produced  signals  with  spectzums  that  peaked  at  frequen¬ 
cies  well  within  the  range  of  the  tape  recorder.  Therefore,  sub¬ 
sequent  experiments  used  Plexiglas  as  the  half-space  and  Styrofoam 
as  the  low-velocity  surface  layer.  It  was  also  anticipated,  how¬ 
ever,  that  because  of  high  attenuation  in  Plexiglas  the  function 
generator  would  not  produce  detectable  signals,  especially  when 
capacitor  pickups  were  being  used.  Therefore,  explosions  were 
used  as  the  source. 

A  layer  of  Styrofoam  10  cm  deep  was  cemented  to  a  Plexiglas 
sheet.  The  pertinent  properties,  plate  velocity  and  density  of 
the  materials  are: 

Vp ( km/sec )  P  ( gm/cm^ ) 

Plexiglas  2.3  1.22 

Styrofoam  0.9  0.03 

A  comparison  of  the  spectrum  of  the  P-wave  as  derived  by 
stated  digital  methods  and  by  use  of  the  Kay  spectograph  is 
shown  in  Figure  17.  The  digital  spectrum  corresponds  to  a  time 
from  onset  to  250  y  sec.  The  analogue  spectrum  was  derived  from 
the  Sonograph  sections  given  in  Figure  12,  The  first  section  in 
Figure  12  represents  the  noise  preceding  the  onset.  The  next 
three  sections  are  separated  at  intervals  of  130  u  secs  (in 
real  model  time),  starting  as  close  to  the  event  as  the  system 
could  be  set.  The  three  sections,  then,  cover  the  same  time 
window  as  the  digital  spectrum.  The  noise  of  section  1  was 
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subtracted  graphically  from  the  other  three.  Then  the  three  sec¬ 
tions  were  added  to  give  an  equivalent  integrated  spectrum.  The 
curves  shown  in  Figure  17  demonstrate  that  the  two  methods  give 
results  ttfhlch  are  quite  close. 

In  preparation  for  the  first  thorough  test  of  the  methods, 
the  transmission  coefficient  through  a  single  low-velocity  layer 
over  a  half-space  was  determined.  The  spectrum  of  the  wave  inci¬ 
dent  at  the  base  of  the  layer  and  that  observed  at  the  top  were 
divided  to  yield  the  frequency-dependent  transmission  coefficient: 


where  I  (u)  and  0  (u)  are  the  observed  input  and  output  spectrums 
respectively.  I  (u)  is  the  digital  curve  in  Figure  17.  The  ex¬ 
perimental  result  is  shown  in  Figure  18.  Three  pronounced  peaks 
occur  at  13,  18,  and  24  kHz. 

Simple  theory  for  the  problem  studied,  e.g.  Backus  ( 1959) # 
predicts  peaks  at  2250,  6750,  11,250,  16,750,  etc.  A  slight 
change  in  the  plate  velocity  in  the  layer,  from  0.9  to  1.1 
km/sec  puts  these  peaks  at  2750,  8250,  13750,  19,250,  24,750, 
etc.  The  actual  speed  in  the  sample  of  Styrofoam  used  will  be 
remeasured  to  determine  if  the  slight  shift  of  these  peaks  can 
be  explained  this  way. 

The  reason  for  the  lack  of  transmission  peaks  above  24  kHz 
is  not  clear.  The  input  spectrum  contains  significant  energy 
out  to  about  50  kHz,  and  yet  these  high  frequencies  are  lost 
in  transmission.  Further  examination  of  the  absorptive  proper¬ 
ties  of  the  Styrofoam  may  clear  this  up. 

» _  -  —  — 


experimental  transmission 
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The  next  step  will  be  to  use  strictly  analogue  techniques 
to  determine  the  transmission  coefficient  and  compare  the  re¬ 
sults.  If  they  are  successful,  the  dipping  bed  problem  will  be 
attacked. 
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B.  Summaries  of  Work  Completed. 

1.  "Love  and  Rayleigh  Wave  Phase  Velocities  over  U.S.  Con¬ 
tinental  Paths,"  by  Michael  L.  Goodwin.  Summary  of  Ph.D. 
dissertation,  September,  1968. 

Surface  waves  of  nine  earthquakes,  with  propagation 
paths  crossing  the  United  States,  were  chosen  for  analy¬ 
sis.  The  Love  and  Rayleigh  waves  were  digitized  from 
seismograms  of  both  WWSSN  and  LRSM  stations.  The  use  of 
both  networks  of  stations  provided  a  greater  number  of 
possible  profiles  in  the  area  approximately  bounded  by 
the  Rocky  Mountains  on  the  west,  and  the  Appalachian 
Mountains  on  the  east.  A  total  of  94  profiles,  shown 
in  Figure  19  were  used  in  this  study. 

The  digitized  data, smoothed,  detrended  and  mean 
removed,  were  Fourier  analyzed  in  the  period  range  be¬ 
tween  10  and  70  seconds,  and  corrected  for  instrumental 
effects.  The  phase  velocities  between  pairs  of  sta¬ 
tions  were  computed,  using  the  difference  in  arrival 
times  of  the  Fourier  components  at  each  of  the  two  sta¬ 
tions  of  the  pair.  An  iteration-least  squares  inversion 
process  was  applied  to  the  observed  Love  wave  and  Ray¬ 
leigh  wave  phase  velocities  to  obtain  the  shear  veloc¬ 
ity  distribution  down  to  about  200  km  deep.  An  earth 
model  of  15  layers  (for  the  200  km  depth)  was  used, 
allowing  more  detailed  results  than  have  been  obtained 
from  thicker  layered  models.  The  calculations  for  the 
Fourier  analysis  and  phase  velocity  and  the  inversion 
procedures  were  run  on  the  IBM  360/50  computer  at 
Washington  University,  St,  Louis,  Missouri. 

The  results  of  the  inversions  show  a  crust  slightly 
greater  than  40  km  in  the  north-central  United  States. 

The  crust  thickens  to  nearer  50  km  under  the  southern 
Canadian  Shield,  Lake  Michigan  and  the  northern  Appala¬ 
chian  Mountains.  The  western  edge  of  the  Great  Plains 
area  has  a  crust  which  is  close  to  50  km  thick,  perhaps 
greater,  and  which  thins  toward  40  km  in  northern  New 
Mexico.  In  southern  Arizona  and  southern  New  Mexico  the 
crust  is  about  40  km  thick,  as  it  is  in  the  central 
United  States  and  Gulf  Coast  areas. 

Except  for  the  Gulf  Coast  area,  most  of  the  pro¬ 
files  indicate  a  low-velocity  channel  in  the  crust  and 
in  the  mantle.  The  low-velocity  channel  in  the  crust 
begins  at  about  21  km  and  extends  to  about  26  km  depth. 
The  low-velocity  channel  in  the  mantle  usually  lies  be¬ 
tween  80  and  180  km  deep.  In  the  places  in  which  the 
low-velocity  channel  is  thinner,  it  often  begins  be'ween 
60  and  80  km  deep.  See  Figures  20  and  21. 

In  the  Gulf  Coast  area,  the  evidence  for  the  low- 
velocity  channels  in  the  crust  and  mantle  indicate  that 
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Figure  20.  Shear  Velocity  vs.  Depth  Profiles  for 

(a)  north-central  United  States  and 

(b)  soutnern  Canadian  Shield -Great 
Lake  Ai’ea. 
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Figure  21.  Shear  Velocity  vs.  Depth  Profile 
for  the  Arizona-New  Mexico  to 
north-central  United  States  Area. 
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the  low-velocity  channel  in  the  crust  disappearsj  and 
the  low-velocity  channel  in  the  mantle  either  disappears 
or  its  depth  is  beyond  the  range  of  the  data.  These  re¬ 
sults  are  summarized  in  the  table. 


TABLE 

SUMMARY  OF  SHEAR  VELOCITY  MODELS* 


PROFILE* 

CRUST 

LVC  IN 

CRUST 

LVC  IN  MANTLE 

TABLE  NO. 

THICK 

•  DEPTH 

VEL. 

DEPTH 

VEL. 

(KM) 

(KM) 

(KM/SEC) 

(KM) 

(KM/ SEC) 

north-central 

u  .  s  •  * 

SOUTHERN 

CANADIAN 

SHIELD* 

GREAT  LAKES 

l  AREAS. 

ATL-HHND* 

8 

40 

26-40 

3.62 

120+ 

4.37 

ATL-GIMA* 

97 

38 

26-38 

3.56 

80+ 

4.48 

ATL-MNN. 

9 

26 

21-26 

3.54 

60-80 

4.40 

atl-wnsd* 

95 

26 

NONE 

- 

100+ 

4.52 

8LA-B0Z* 

92 

38 

21-26 

3.49 

.  100-160 

4.28 

BLA-WNSD* 

94 

38 

21-26 

'  3.56 

* 

BRPA-BOZ* 

77 

32 

21-26 

3.42 

60+ 

4.28 

FRMA-BLA* 

20 

26 

NONE 

- 

80-160 

4.30 

GIMA-8LA* 

23 

32 

NONE 

- 

60-140 

4.10 

GIMA-GEO. 

60 

32 

NONE 

- 

80+ 

4.18 

SCP-GIMA* 

15 

26 

21-26 

3.51 

100+ 

4.41 

GIMA-SCP* 

57 

38 

21-26 

3.57 

100+ 

4.36 

RYND-OGD* 

56 

26 

21-26 

3.41 

100-160 

4.22 

BRPA-MNN  * 

87 

50 

NONE 

- 

* 

- 

ARIZONA-NEW  MEXICO 

TO 

NORTH-' 

CENTRAL 

U.S.  AREA. 

TDNM-MDS* 

41 

40 

21-26 

3.62 

60-100 

4.11 

TDNM-MNN * 

49 

38 

NONE 

- 

100-120 

4.37 

TUC-AAM* 

44 

38 

NONE 

- 

80-160 

4.13 

TUC-MDS, 

51 

32 

21-26 

3.56 

100+ 

4.16 

TUC-NNN* 

25 

38 

NONE 

- 

* 

- 

TUC-EBMT* 

45 

4o 

NONE 

- 

* 

DRCO-EBMT  * 

48 

4q- 

50 

NONE 

- 

12o+ 

4.4o 

TABLE  (CONTINUED) 


profile* 

CRUST 

LVC  IN 

CRUST 

LVC  IN  MANTLE 

table  no. 

THICK, 

DEPTH 

VEL. 

DEPTH 

VEL. 

(KM  J 

(KM) 

(KM/SEC) 

'  (KM) 

(KM/SEC 

PENNSYLVANIA  TO  ARIZONA-NEW  MEXICO  AREA, 

OGD-ALQ  » 

63 

26 

NONE 

.  *"» 

100-120 

4.37 

SCP-ALQ. 

73 

40-50 

NONE 

120-140 

4.46 

SCP-TUC. 

76 

32 

NONE 

'  «. 

100+ 

4.52 

TUC-BLA. 

43 

40 

NONE 

- 

180+ 

4.22 

TDNM-SCP* 

50 

40 

NONE 

- 

80+ 

4.14 

SKTX-BLWV, 

38 

38 

NONE 

160+ 

4.61 

GEO-TUC. 

74 

32 

NONE 

- 

100-120 

4.34 

SKTX-DHNY, 

30 

38 

21-26 

3,51 

60+ 

4.76 

SKTX-OGD* 

36 

32 

21-26 

,3.40 

..  60-80 

4.59 

FLO-TUC* . 

62 

40 

NONE 

- 

100+ 

4.12 

flo-alq. 

64 

38 

21-26 

3.51 

100+ 

4.12 

SOUTH-CENTRAL 

U.S.  AREA, 

.  .. 

ATL-DRCO* 

80 

32-38 

NONE 

_ 

80-180 

4.13 

ATL-GOL* 

83 

32-38 

21-26 

3.54 

80-180 

4.20 

GOL-ATL, 

38 

32 

21-26 

3.52 

60+ 

4.22 

ATL-DUG  » 

80 

26 

21-26 

3.43 

60+ 

4.10 

OXF-RTNM* 

84 

38 

21-26 

3.53 

60+ 

4.02 

SHA-ALQ* 

86 

38 

21-26 

3',  54 

60-120 

4.13 

TUC-OXF. 

29 

38 

NONE 

- 

140+ 

4.58 

TUC-EUAL* 

31 

40 

NONE 

- 

SEE  TEXT 

• 

TUC-SHA* 

33 

33 

21-26  r 

3.56 

100+ 

4.39 

TUC-OAL* 

35 

38 

NONE 

- 

60-160 

4.21 

SHA-LUB* 

91 

33 

21-26 

3.50 

60+ 

4.04 

LUQ-ATL* 

42 

38 

21-26 

3.49 

60+ 

4.00 

LUB-OXF, 

54 

38 

21-26 

3,55 

60-160 

4.29 

ATL-OXF. 

79 

38 

NONE' 

- 

100+ 

4.51 
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"Decoupling  and  Source  Function  for  Explosions  in  a 
Three-Dimensional  Model,"  by  Dhari  Saaid  Bahjat,  Digest 
of  Ph.D.  dissertation.  May,  1968. 

A  three-dimensional  model  made  of  Plaster  of  Paris 
has  been  used  to  investigate  some  conditions  of  an  ex¬ 
plosion  which  control  the  seismic  signal  generated  by 
such  a  source.  This  study  is  based  on  the  idea  that 
the  seismic  signal  generated  by  an  explosion  can  be  re¬ 
duced  by  detonating  the  explosive  in  a  cavity.  Spherical 
charges  of  chemical  explosive  are  used  to  generate  the 
seismic  signal.  The  seismic  waves  are  recorded  at  the 
surface  of  the  model  by  calibrated  capacitor  receivers. 
The  compressional  wave  velocity  and  absorption  in  Plaster 
of  Paris  are  determined. 

The  factors  which  affect  the  seismic  signal  gen¬ 
erated  by  the  tamped  explosion  of  a  particular  explo¬ 
sive  material  are  the  yield  of  the  explosion,  the 
depth  of  the  explosion  and  the  medium  in  which  the  ex¬ 
plosion  is  buried.  Data  are  collected  for  nine  charge 
masses  in  the  range  from  20  to  162  milligrams.  The  re¬ 
sults  show  that  the  peak  amplitude  of  the  compressional 
wave  depends  on  0.72  power  of  the  yield.  Also  it  is 
observed  that  the  pulse  width  varies  with  the  size  of 
the  charge.  The  dependence  is  0.6  power  of  the  yield. 

Two  series  of  experiments,  one  using  55  milligrams 
charges,  the  other  100  milligrams  charges,  have  been 
conducted  to  investigate  the  variation  of  the  peak  am¬ 
plitude  of  compressional  waves  with  the  cavity  radius. 

The  cavity  radii  are  varied  from  0.31  to  3.2  centi¬ 
meters  . 

The  amplitude  of  the  compressional  wave  for  a 
tamped  shot  is  smaller  than  the  one  for  cavity  shots 
for  the  range  of  cavity  sizes  investigated.  In  the 
case  of  tamped  shots,  most  of  the  energy  released  by 
the  explosion  is  consumed  in  crushing  the  surrounding 
medium,  and  as  waste  heat  deposited  in  the  source  re¬ 
gion.  The  amplitude  of  the  compressional  wave  in¬ 
creases  as  the  cavity  radius  Increases  as  long  as  the 
pressure  applied  to  the  surface  of  the  cavity  exceeds 
the  compressive  strength  of  the  medium.  The  maximum 
amplitude  is  obtained  for  a  cavity  radius  such  that 
the  stress  transmitted  through  the  wall  of  the  cavity 
to  the  medium  is  equal  to  the  compressive  strength  of 
the  medium.  The  radius  of  the  cavity  in  which  the 
pressure  applied  on  the  surface  of  the  cavity  is  equal 
to  the  compressive  strength  of  the  medium  is  a  func¬ 
tion  of  the  charge  mass.  The  amplitude  decreases  after 
reaching  a  maximum  because  of  the  reduction  in  the  ef¬ 
fective  source  region  radius.  Therefore,  the  factor 
by  which  the  amplitude  is  reduced  for  a  given  increase 
in  cavity  radius  is  independent  of  the  yield  of  the 
explosion. 


83 


The  energy  calculations  verify  that  only  a  small 
fraction  of  the  released  energy  Is  radiated  Into  seis¬ 
mic  waves.  The  seismic  energy  radiated  In  the  P-wave 
is  computed  from  the  early  part  of  the  seismograms. 

It  follows  the  same  general  pattern  for  Increase  and 
then  decrease  with  Increasing  radius  that  the  peak  am¬ 
plitude  does.  The  radiated  seismic  energy  for  a  tamped 
shot  Is  shown  to  be  less  than  for  the  cavity  shots  for 
this  range  of  cavity  sizes. 

The  Idea  that  for  some  cavities  the  radius  of  the 
source  region  Is  larger  than  the  radius  of  the  given 
cavity  Is  checked  by  observing  the  change  In  the  pulse 
width  of  the  first  arrival.  A  minimum  pulse  width 
occurs  at  larger  cavity  radii  for  bigger  charges.  This 
Is  explained  very  well  by  the  linear  theory  of  a  cavity 
In  an  Infinite  medium. 

The  maximum  stress  on  the  surface  of  the  equiva¬ 
lent  cavity  Is  determined  In  two  Independent  ways. 
First,  starting  from  the  point  of  explosion  and  pro¬ 
ceeding  outward  from  the  source  region  to  the  elastic 
region,  and  secondly,  starting  from  seismic  waves  ob¬ 
served  at  the  surface  of  the  model  and  moving  Inward 
to  the  elastic  boundary  around  the  source.  The  values 
of  the  peak  stress  obtained  from  equalization  of  re¬ 
corded  seismic  waves  to  the  surface  of  the  cavity  are 
In  satisfactory  agreement  with  the  computed  Input 
stress  to  the  medium  for  cavities  which  respond  elas¬ 
tically.  For  overdriven  cavities  the  Input  stress 
cannot  be  compared  with  that  of  calculated  one  from 
the  seismic  zone  unless  the  decay  of  the  stress  In 
the  fracture  region  Is  known. 

"Exact  Solution  for  the  Displacement  Components  In 
P1P2P1*  (sls2sl)sv»  (S1S2S1^SH  Head  Waves  Due  to  an 
Impulsive  Double  Couple  Source,"  by  Umesh  Chandra, 

In  press.  Bull.  Seism,  Soc.  Am,,  39  (1969) 

Abstract 

Exact  Integral  expressions  for  the  displacement 
components  as  a  function  of  azimuth.  In  P]^^* 

(S1SgS1)SY  t  and  (S^SgS^JgH  head  waves  due  to  an  Im¬ 
pulsive  double  couple  source  of  arbitrary  orientation 
In  two  semi-infinite  media  separated  by  a  plane  hori¬ 
zontal  boundary,  and  the  effect  of  free  surface  on  It, 
have  been  obtained.  The  time  variation  of  the  source 
Is  assumed  to  be  represented  by  the  Heaviside  unit 
function.  The  response  due  to  an  a’^bltrary  source 
excitation  may,  however,  be  obtained  from  this  solu¬ 
tion  by  superposition.  Necessary  formulae  for  the 
operational  interpretation  of  Integral  equations 
occurring  In  the  problem  are  given. 


Illustrative  examples  showing  the  \-javeform  in 

P1P2P1'  (sls2Sl)sV*  and  (S1S2S1^SH  head  wave  dlsPlace“ 
ment  components  for  a  particular  fault  and  a  suitably 

chosen  earth  model  are  given.  Plots  of  the  projection 
of  particle  trajectory  for  P  waves  in  ur-u2  plane,  and 
for  S  waves  in  u^-ur  and  u^-Ug  planes  are  given  for 
different  azimuths;  where  ur,  u  ,p  ,  uz  are  the  compo¬ 
nents  of  displacement  in  a  cylindrical  coordinate  sys¬ 
tem  (r, f  ,  z) . 
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